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Chapter  1 


Introduction 
1.1  Motivation 

Engineering  models  today  are  mainly  developed  within  a  3D  computer  aided 
design  (CAD)  software,  for  the  reasons  that  the  CAD  models  are  unambiguous, 
visualization  and  manufacturing-planning  is  easier  within  the  environment,  and  3D 
model  data  exchange  is  well  established.  The  3D  models  created  using  these  CAD 
systems  often  contain  an  abundant  amount  of  features.  Fig.  1.1  illustrates  an  air¬ 
craft  engine  with  thousands  of  features.  These  features  are  necessary  for  a  wide 
variety  of  life  cycle  needs  including  safety,  functionality,  aesthetics,  and  manufac¬ 
turability. 

The  generated  3D  CAD  models  are  used  for  downstream  applications  within 
the  design  and  manufacturing  process.  Finite  element  analysis  (FEA),  is  one  of  the 
downstream  applications,  has  emerged  as  an  important  engineering  analysis  tool 
and  is  extensively  used  for  structural,  thermal,  and  modal  analyses  to  ensure  that 
the  design  of  the  model  will  meet  the  specifications  required  from  the  operation 
of  the  model.  FEA  is  a  numerical  technique  for  calculating  approximate  solutions 
of  complicated  equations  used  to  calculate  the  analysis  results.  It  does  this  by 
generating  a  mesh,  grid  of  nodes,  over  the  entire  model.  The  density  of  the  grip 
varies  throughout  the  model,  using  having  a  higher  density  at  hllets,  corners,  and 
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high  stress  areas.  The  mesh  contains  the  structural  and  material  properties,  which 
define  how  the  structure  reacts  to  the  conditions  implemented  within  the  analysis. 
At  each  node  an  equilibrium  equation  is  generated  and  the  equations  are  solved  to 
obtain  the  results  in  the  analysis. 

Importing  the  detailed  3D  models,  such  as  the  one  shown  in  Fig.  1.1,  into 
the  FEA  software  often  presents  several  challenges.  First,  the  meshing  time  rapidly 
increases  with  the  complexity  of  the  model.  Second,  it  may  take  a  long  time  to 
analyze  the  large  FEA  mesh  that  results  from  the  detailed  CAD  models.  Third,  the 
presence  of  a  small  feature  may  lead  to  a  poor  quality  mesh  and  hence  may  lead 
to  erroneous  analysis  results.  In  extreme  cases,  the  automated  meshing  may  fail 
on  very  complex  models,  requiring  manual  meshing  of  the  model.  Utilizing  a  more 
powerful  computer  can  not  solve  all  the  issues  associated  with  analyzing  complex 
CAD  models.  Meshless  methods  have  been  proposed  as  an  alternative  to  FEA, 
which  use  the  geometry  of  the  object  directly  for  the  calculations.  These  methods 
are  not  as  widely  used  within  the  community. 
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Figure  1.1:  Napier  Sabre  24  cylinder  aicraft  engine  [1] 
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To  circumvent  the  above  described  problems,  experienced  engineering  analysts 
rarely  perform  meshing  on  the  original  complex  models  containing  all  the  design 
features.  In  most  situations,  the  model  is  simplihed  before  proceeding  with  meshing. 
Simplihcation  carefully  examines  the  model  and  suppresses  the  features/parts  that 
are  not  likely  to  influence  the  analysis  results,  but  may  signihcantly  improve  the 
computational  performance  both  in  terms  of  mesh  size  and  mesh  quality.  Currently, 
the  simplihcation  of  models  is  mainly  performed  manually,  which  can  take  a  long 
time  to  perform  on  a  complex  model  and  can  be  tedious  in  nature.  Also,  most  of  the 
state-of-the-art  mesh  generation  algorithms  fail  to  generate  meshes  for  those  features 
that  are  smaller  than  the  FE  size.  In  these  situations,  CAD  model  simplihcation  is 
essential  to  producing  FE  models  adaptable  to  analysis.  Certainly,  automation  in 
this  task  can  improve  analysts’  productivity.  Hence,  automated  feature  suppression 
has  been  acknowledged  as  a  useful  tool  in  CAD-FEA  integration  [26]. 

Automated  optimization  of  a  CAD  model  for  analysis  can  be  a  highly  benehcial 
tool  within  the  industry.  As  mentioned  previously,  the  models  are  currently  simpli¬ 
hed  manually  for  hnite  element  analysis.  Automating  this  process  can  potentially 
decrease  the  amount  of  time  required  for  simplihcation,  suppress  all  non-critical  fea¬ 
tures/parts,  and  retain  all  critical  fetaures/parts;  ultimately  speeding  up  the  process 
and  reducing  the  amount  of  errors  involved  within  the  process.  For  automation  of 
this  simplihcation  process,  numerous  factors  will  have  to  be  accounted  for  to  cor¬ 
rectly  classify  a  feature/part  as  critical  or  non-critical  to  the  analysis.  A  few  of  these 
factures  include  dimensions,  material,  and  purpose. 

Two  diherent  types  of  approaches  can  be  taken  to  perform  automated  model 
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simplification  for  FEA.  The  first  approach  is  geometry  based  and  utilizes  the  geo¬ 
metric  operators  to  simplify  the  parts  [40].  Unfortunately,  encoding  the  preferences 
of  an  analyst  in  this  approach  is  very  difficult  and  can  lead  to  erroneous  and  in¬ 
consistent  results  [30].  The  second  approach  is  knowledge  based.  In  this  case,  the 
preferences  and  knowledge  used  by  the  analyst  are  explicitly  encoded  as  simplihca- 
tion  rules.  This  approach  captures  analysts’  knowledge  in  terms  of  rules  that  express 
the  characteristics  of  features/parts  that  should  be  considered  for  suppression. 

1.2  Thesis  Goals  and  Scope 

The  goals  and  scope  of  this  thesis  revolve  around  developing  a  simplihcation 
tool  for  complex  CAD  models  in  preparation  for  FEA.  The  FEA  types  that  are 
focused  on  within  the  framework  are  structural  and  modal.  This  tool  is  divided 
into  two  primary  parts:  part  model  simplihcation  through  feature  suppression  and 
assembly  model  simplihcation  through  part  suppression  and  part  simplihcation.  For 
development  of  these  tools  a  knowledge  based  approach  has  been  utilized. 

1.2.1  Part  Simplification 

Part  model  simplihcation  consists  of  suppressing  non-critical  features  within 
a  CAD  model.  The  determination  of  the  non-critical  (suppressible)  features  relies 
on  the  diherent  attributes  of  not  only  the  features  themselves,  but  also  of  the  entire 
part  model  and  analysis.  Some  of  these  attributes  include  the  feature  type,  feature 
dimensions,  feature  distances  to  boundary  conditions,  analysis  type,  and  part  di- 
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mensions.  The  dimensions  of  the  part  is  required,  so  that  the  relative  size  of  the 
feature  to  the  overall  part  can  be  calculated.  This  is  crucial,  since  part  sizes  are 
not  limited  and  a  feature  that  is  seen  as  dimensionally  small  within  one  part  may 
be  quite  large  within  another.  Fig.  1.2  illustrates  an  example  of  a  CAD  model  with 
features  capable  of  being  suppressed  and  those  features  that  cannot  be  suppressed. 


(a)  Boundary  conditions 


(b)  Suppressibility  of  features 


Figure  1.2:  Suppressibility  of  a  feature:  The  decision  to  suppress  or  preserve  a  feature 
strongly  depends  upon  the  analysis  context  and  the  application.  For  a  given  loading 
and  boundary  conditions,  “Insuppressible  features”  would  affect  the  FEA  analysis 
whereas  the  “suppressible  features”  would  not. 


1.2.2  Assembly  Simplification 

Assembly  model  simplihcation  consists  of  suppressing  or  simplifying  part  mod¬ 
els  within  a  CAD  model.  Similarly  to  that  of  the  part  approach,  parts  can  be  labeled 


5 


as  either  non-critical  (suppressible)  or  critical  (non-suppressible);  although  unlike 
the  part  appproach,  there  can  also  be  semi-critical  parts  within  the  assembly  model. 
Semi-critical  are  the  parts  that  must  remain  in  the  model  for  analysis,  but  the  over¬ 
all  complexity  of  the  part  is  not  required;  this  can  either  consist  of  suppressing  the 
non-critical  features  or  simplifying  the  overall  geometry  of  the  part. 

1.3  Organization 

Chapter  2  is  a  survey  of  the  current  literature  within  the  held  of  model 
simplihcation,  which  is  not  only  limited  to  the  simplihcation  techniques  specihed 
for  FEA  model  preparation.  This  chapter  is  arranegd  into  four  collective  sections; 
model  simplihcation  based  upon  surface  operators,  volume  operators,  dimension 
reduction  operators,  andexplicit  feature  operators.  Each  section  is  then  further 
split  up  into  the  specihc  methods  types  of  the  simplihcation. 

Chapter  3  presents  a  model  simplihcation  framework  for  CAD  part  models, 
by  suppressing  non-critical  features.  In  this  framework  a  decision-tree  based  rep¬ 
resentation  has  been  developed  for  implementing  the  suppression  rules  on  features; 
such  as  holes,  rounds,  and  chamfers.  Eeature  suppression  demonstrations  from  an 
expert  analysts  provide  a  useful  way  of  acquiring  application  dependent  feature 
suppression  rules.  This  method  does  not  burden  the  human  expert  with  explicitly 
cataloging  feature  suppression  rules.  Instead,  feature  suppression  rules  are  extracted 
through  demonstrations,  in  an  unobtrusive  fashion.  These  application  dependent 
rules  can  be  learned  by  using  decision  tree  learning  algorithms.  We  have  developed 
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utilities  in  the  Pro/Engineer  CAD  system  to  identify  features  that  meet  the  learned 
rules  and  suppress  them. 

Chapter  4  presents  a  model  simplihcation  framework  for  CAD  assembly  mod¬ 
els,  by  suppressing/simplifying  non-critical  parts.  In  this  framework  utilities  have 
been  developed  within  the  Pro/Engineer  CAD  system  to  identify  and  suppress  parts 
with  a  specified  material  or  PART_NAME  value. 

Chapter  5  discusses  the  intellectual  contributions  of  the  research  and  possible 
areas  of  future  research  to  further  extend  the  frameworks  discussed  within  the  thesis. 
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Chapter  2 


Related  Work 

2.1  Overview 

Model  simplification  has  been  a  focus  within  the  modeling  area  for  a  while 
now,  consequently  there  is  a  rich  body  of  literature  within  this  held.  As  not  all 
of  the  literature  can  be  reviewed  within  this  thesis,  I  will  focus  on  some  of  the 
more  prevalent  proposed  techniques.  For  organization  purposes  the  section  follows 
the  classihcation  within  a  survey  paper  by  Thakuer  et  ah  [40],  specihcally:  surface 
based  operators,  volume  based  operators,  dimension  reduction  based  operators,  and 
explicit  feature  based  operators. 

2.2  Techniques  Based  on  Surface  Entities 

This  section  reviews  the  techniques  for  which  the  simplihcation  is  based  upon 
the  surface  entities  and  are  organized  into  three  classihcations;  boundary  loop  de¬ 
composition,  low  pass  hltering,  and  face  cluster  based  simplihcation. 

2.2.1  Boundary  Loop  Decomposition 

Sun  et  al.  [39]  presents  a  simplihcation  method  for  B-rep  models  using  region 
suppression.  First,  the  simplihed  region  (SR)  and  its  boundary  loop  must  be  deter- 


mined  and  classified.  A  hybrid  method,  both  automated  and  manual,  was  adopted 
for  determination  of  the  SR  and  the  faces  which  are  to  simplified  are  identified  us¬ 
ing  their  area  ratio.  After  identifying  the  SR,  the  boundary  loops  are  formed  from 
the  edges  of  the  region.  Each  SR  is  classified  according  to  whether  it  has  single 
or  multiple  boundary  loops.  Multiple  loop  boundaries  must  be  reduced  to  a  single 
loop,  either  by  dealing  with  the  loops  separately  or  constructing  a  bridge  between 
the  two  separate  loops  to  combine  the  loops  into  one. 

A  divide  and  conquer  strategy  is  then  used  to  decompose  the  loops.  This 
recursively  divides  the  boundary  loops  into  multiple  sub-loops  by  connecting  two 
non-adjacent  vertices  with  an  edge,  if  there  are  either  critical  edges  or  boundary 
edges  lying  on  the  same  underlying  curve.  The  iteration  process  continues  until 
the  no  co-defined  edge  pairs  can  be  identified.  The  sub-loops  are  then  further 
decomposed  into  co-surface  loops. 

Upon  the  decomposition  of  the  boundary  loops,  the  SRs  are  deleted  from  the 
model.  The  hole  is  fixed  by  REPAIR  and  PATCH  operators.  The  REPAIR  operator 
modifies  the  neighboring  faces  boundary  and  the  PATCH  operator  constructs  a  new 
face  and  stitches  the  faces  to  the  model.  Both  operations  are  possible  through  the 
use  of  a  co-surface  loop.  To  validate  the  B-rep  model  duplicate  edges,  vertices, 
and  faces  are  merged.  Advantages  of  the  technique  described  is  that  it  does  not 
require  Boolean  operators  and  the  REPAIR  operator  is  insensitive  to  different  types 
of  features. 
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2.2.2  Low  Pass  Filtering 


Lee  and  Lee  et  al.  [29]  proposed  a  semi-automatic  low  pass  filtering  technique 
using  Discrete  Fourier  Transform.  For  implementation  of  the  method  the  CAD 
model  must  hrst  be  digitized,  converted  into  a  2D  digital  image  of  resolution  512 
X  512.  A  surface  can  then  be  plotted  from  the  pixels,  having  a  value  of  -1.0,  0,  or 
1.0.  Upon  application  of  the  Fourier  transform  to  the  surface  function,  the  surface  is 
expressed  in  the  frequency  domain.  Within  the  frequency  domain  the  low  frequency 
terms  represent  the  overall  shape,  whereas  the  high  frequency  terms  represent  the 
detailed  features.  Removal  of  the  high  frequency  components  generates  a  low  pass 
hltered  (LPF)  model  of  the  original  model.  The  adverse  result  of  the  LPF  model  is 
that  all  sharp  corners  were  rounded.  Removal  of  the  sharp  edges  is  unfavorable  in  the 
analysis  model,  since  sharp  edges  tend  to  acquire  a  higher  stress  concentration.  To 
retain  the  sharp  edges  within  the  simplihed  model,  the  boundary  elements  average 
distance  between  the  original  model  and  LPF  model  is  evaluated  as  a  metric  of 
complexity.  If  the  metric  is  below  a  given  threshold  the  entity  is  considered  to  be 
detailed  and  if  the  metric  is  below  the  given  threshold  the  entity  is  considered  to  be 
to  belong  to  the  overall  shape.  The  method  is  only  developed  for  2D  images  and  if 
extended  to  a  3D  image  it  is  limited  in  terms  reconstructing  the  removed  faces. 

2.2.3  Face  Cluster  Based  Simplification 

Inoue  et  al.  [18]  proposed  a  face  clustering  technique  for  CAD  model  simplih- 
cation.  The  hrst  step  within  the  process  is  to  decompose  the  model  into  faces.  From 
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this  decomposition  a  planar  adjacency  graph  can  be  constructed;  the  graph  nodes 
represent  the  faces,  which  are  connected  by  arcs  representing  the  adjacency  of  the 
face  pairs.  To  determine  the  nodes  to  be  clustered,  a  weight  for  each  arc  is  calculated 
based  on  several  geometric  indices.  The  geometric  indices  attempt  to  take  the  area 
size,  boundary  smoothness,  and  region  flatness  into  consideration.  The  equations 
for  these  indices  are  presented  within  the  paper  for  reference. 

Upon  calculation  of  the  arc  weights,  the  arcs  are  prioritized  for  node  clustering 
by  ranking  the  weight  values  highest  to  lowest.  After  prioritizing  the  arcs,  a  merge- 
ability  test  is  completed  using  the  geometric  indices.  The  top  weighted  arc  that 
passes  the  mergeability  test  is  then  contracted,  merging  the  two  connected  nodes.  If 
the  two  nodes  both  had  an  connected  with  a  similar  node,  the  two  arcs  are  combined 
into  one.  This  process  is  then  completed  for  the  entire  adjacency  graph,  simplifying 
the  original  model.  Fig.  2.1  display  the  process  of  removing  the  arcs  and  merging 
the  nodes.  For  the  example  discussed  within  the  paper,  a  model  containing  1742 
individual  faces  was  simplihed  to  a  hnal  model  containing  only  34  faces. 


Figure  2.1:  Face  clustering  [18] 


Sheffer  et  ah  [2]  expanded  upon  this  technique  proposed  by  Inuoe.  The  clus¬ 
tering  algorithm  is  maintained  within  the  new  process,  although  different  geometric 
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indices  are  calculated.  The  new  geometric  indices  quantify  the  boundary  preser¬ 
vation,  region  size,  region  smoothness,  and  simplicity  of  the  boundary  shape.  The 
main  advantage  of  these  indices  is  that  they  can  also  account  for  non-planar  regions. 

In  addition  to  the  new  indices,  Sheffer  introduced  a  collapsibility  test  to  de¬ 
termine  faces  that  can  be  symetrically  divided  between  adjacent  clusters,  after  the 
clustering  process.  This  division  is  completed  using  virtual  topology  operators  to 
merge  and  collapse  the  faces.  In  this  process  the  face  cluster  found  to  be  divisible 
from  the  collapsibility  test  is  then  split  into  a  number  of  faces  equal  to  the  number 
of  adjacent  neighboring  faces.  The  new  faces  are  then  merged  to  the  respective 
neighboring  face  from  which  an  edge  is  shared.  This  further  reduces  the  number 
faces  within  the  model. 

In  a  similar  technique  to  the  cluster  division  by  virtual  topology  operators, 
Dey  et  ah  [10]  introduced  a  framework  to  eliminate  the  feature  effects  at  mesh 
level.  This  technique  identihes  mesh  entities  with  per  quality  metrics,  based  upon 
the  dihedral  angle  measurement  and  aspect  ratio  of  the  mesh.  After  identihcation, 
the  suppression  of  these  entities  is  ensured  to  maintain  geometric  validity  of  the 
modihed  mesh  and  the  entities  are  suppressed  accordingly. 

2.3  Techniques  Based  on  Volumetric  Entities 

This  section  reviews  the  techniques  for  which  the  simplihcation  is  based  upon 
the  surface  entities  and  are  organized  into  two  classihcations;  voxel  based  simplih¬ 
cation  and  effective  volume  simplihcation. 
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2.3.1  Voxel  Based  Simplification 


Andujar  et  al.  [3]  proposed  a  topology-reducing  simplification  algorithm  called 
Trihedral  Discretized  Polyhedra  Simplification  (TDPS).  This  algorithm  contains 
three  steps:  discretization,  reconstruction,  and  face  reduction.  Maximal  division 
classical  octree  (MDCO)  is  used  to  discretize  the  model.  MDCO  arranges  eight 
octants,  obtained  by  recursively  subdividing  the  cubic  universe  into  an  8-ary  tree. 
The  leaves  of  the  tree  falling  completely  in  or  out  of  the  model  volume  are  labeled  as 
either  black  or  white,  respectively,  and  the  leaves  containing  part  of  the  boundary, 
partially  in  or  out  of  the  model,  are  labeled  as  grey.  The  grey  leaves  are  then  further 
subdivided  until  the  parts  are  either  black  or  white  or  the  tree  reaches  a  certain  level. 

Reconstruction  step  involves  converting  the  MDCO  octree  representation  into 
a  3D  digital  picture,  set  of  points  arranged  in  a  regular  grid,  and  extracting  an 
isosurface  from  the  digital  picture.  To  compute  the  digital  picture  the  MDCO  is 
traversed  to  find  the  BTC  nodes,  boundary  nodes.  Each  BTG  node  is  subdivided 
into  eight  octants  and  the  color  at  each  octant  corner  is  then  labeled  as  either  white 
or  black,  depending  on  the  color  of  the  node  in  contact  with  the  corner.  From  this 
digital  picture  a  surface,  R,  is  extracted  based  upon  an  error  bound  in  terms  of  the 
solid-Hausdorff  distance,  satisfying  the  condition  that  all  white  nodes  are  completely 
out  of  R  and  all  black  nodes  are  completely  inside  the  R.  The  face  reduction  step 
further  reduces  the  geometric  complexity  of  R,  using  an  edge  collapse  technique. 

The  authors  report  that  TDPS  has  several  advantages  over  other  model  simpli¬ 
fication  techniques;  suitable  for  both  parts  and  assemblies,  it  is  not  affected  by  input 
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surface  degeneracies,  and  produces  a  two-manifold  solids  even  with  non-manifold  in¬ 
puts.  Although  there  is  a  limitation  of  TDPS  in  that  it  is  not  intended  for  small 
error  thresholds. 

2.3.2  Effective  Volume  Based  Simplification 

Frameworks  for  multiresolution  representations  of  solid  models  producing  sim- 
plihed  models  at  various  levels  of  detail  (LODs)  have  been  extensively  researched. 
[26,  28,  6,  25]  The  focus  of  the  research  can  be  split  into  two  separate  approaches. 
The  first  approach  uses  conventional  solid  boundary  representations  (B-rep)  [6]  and 
the  second  approach  uses  non-manifold  B-reps  [28].  Conventional  solid  B-reps  is  ad¬ 
vantageous  since  it  can  be  implemented  into  commercial  3D  CAD  systems,  due  to 
the  same  data  structures;  although  the  computational  cost  of  the  boundary  evalua¬ 
tion  is  very  high,  causing  the  method  to  be  time  consuming.  Unlike  the  conventional 
B-reps,  non-manifold  B-reps  are  less  time  consuming  and  therefore  will  be  explored 
further. 

The  non-manifold  model  is  capable  of  representing  any  combination  of  wire¬ 
frame,  surface,  solid,  and  cellular  models  in  a  unihed  data  structure.  It  consists 
of  three  key  modules:  a  feature-based  modeling  module,  a  feature-based  idealiza¬ 
tion  module,  and  a  non- manifold  topological  (NMT)  modeling  kernel.  The  modules 
manage  the  database  of  the  features,  perform  the  detail  removal  and  dimension 
reduction,  and  manipulate  the  geometric  models,  respectively.  The  detail  removal 
within  the  idealization  module  consists  of  three  steps  to  simplify  the  model.  First  all 
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idealization  features  with  a  design  application  domain  are  extracted  from  the  master 
model.  Then  according  to  the  given  LOD  criterion,  most  often  the  feature  volume, 
the  extracted  idealization  features  are  rearranged.  If  the  criterion  is  based  upon 
decreasing  volume,  features  with  a  volume  below  the  given  threshold  are  considered 
for  simplihcation.  Finally,  the  LOD  is  selected  manually  to  obtain  the  appropriate 
simplihed  model.  [26] 

Both  the  conventional  and  non-manifold  B-rep  approaches  have  a  feature  ar¬ 
rangement  issue.  Rearrangement  of  the  features  to  generate  the  various  LODs  of 
the  model  will  possibly  result  in  a  shape  different  from  that  of  the  original  model 
shape,  since  the  additive  and  subtractive  Boolean  operations  are  not  commutative. 
Lee  et  ah  [25]  and  Choi  et  ah  [6]  introduce  a  delta  volume,  redehned  volume,  if 
required,  of  the  subtractive  feature  operations  to  ensure  model  similarity.  The  rear¬ 
rangement  algorithm  applied  within  these  approaches  combines  all  additive  features 
to  generate  the  lowest  resolution  model  and  the  delta  volumes  are  subtracted  from 
the  model  in  descending  order  of  volumes  to  generate  the  higher  resolution  models. 
This  method  is  only  applicable  for  the  specihed  feature  rearrangement,  not  for  any 
arbitrary  rearrangement. 

Lee  [28]  expanded  upon  technique  by  introducing  the  concept  of  the  effective 
volume  of  a  feature,  fractional  or  entire  volume  of  the  feature  to  provide  the  same 
result  in  spite  of  Boolean  arrangement.  This  volume  is  calculated  using  Boolean 
algebra.  Let  Vi  denote  the  volume  of  the  solid  primitive  of  a  feature  Xj  denote 
a  Boolean  operation,  and  denote  the  resulting  model  obtained  by  applying  n 
Boolean  operations  between  n-|-l  solid  models: 
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If  the  Boolean  operation  XjVj  is  moved  to  the  mth  position  can  be 
represented  as  follows: 
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Figs.  2. 2-2. 3  illustrate  the  use  of  Eqns.  2. 1-2. 3  on  a  model  generated  by  ap¬ 
plying  Eve  form  features  [27].  The  original  model  is  shown  in  Fig.  2.2.  If  we  alter 
the  order  of  features  from  Fq  — )■  Fi  — )■  F2  — )■  F3  — )■  F4  to  Fq  — )■  F2  — )■  F3  — )■  F4  — )■  Fi 
by  moving  Fi  to  the  last  position,  the  Boolean  operation  sequence  changes  from 
Vo  —  V1UV2  —  V3UV4  to  V0UV2  —  V3UV4  —  Vi,  where  Vi  is  the  respective  volume  of  F*. 
When  just  reorganizing  the  form  features  Vi  will  just  remove  the  volume  for  which 
it  encompasses,  i.e.  removing  V2  and  the  connection  of  V2  to  Vq.  Although,  if  we 
apply  Eqns.  2. 1-2. 3  the  effective  volumes  of  Fq,  Fi,  F2,  F3,  and  F4  are  calculated  to 
be  Vo,  Vi  —  V2  —  V4,  V2,  V3,  and  V4,  respectively.  Fig.  2.3  displays  the  rearranged  fea¬ 
ture  model  using  the  effective  volumes  of  the  features,  which  is  identical  the  overall 
original  model.  This  technique  can  allow  for  arbitrary  rearrangement  of  the  features 
and  generate  the  same  final  shape  as  the  original  model. 
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FO(Base)  F1(Blind_Hole) 


Figure  2.2:  Original  model  construction  [28] 


LOD  =  4 


Figure  2.3:  Reconstruction  of  the  model  by  rearranging  the  features  and  using  the 


effective  volumes  [28] 


2.4  Techniques  Based  on  Explicit  Features 


This  section  reviews  the  techniques  for  which  the  simplihcation  is  based  upon 


the  explicit  features  and  are  organized  into  three  classihcations;  prismatic  feature 

17 


simplification,  blend  simplification,  arbitrary  shaped  feature  simplification. 

2.4.1  Prismatic  Feature  Simplification 

Dabke  et  al.  [8]  implemented  DESIDE-X,  a  rule-based  software  package  for 
suppressing  features,  which  relies  upon  a  constructive  solid  geometry  (CSG)  tree 
and  is  specihcally  for  modal  analysis.  Through  the  use  of  a  feature  by  feature 
looping  mechanism,  each  feature  in  the  tree  is  analyzed  and  determined  to  be  either 
signihcant  or  non-signihcant.  The  signihcance  of  the  feature  is  evaluated  based  upon 
a  rule  base  that  takes  into  account  different  criteria;  such  as,  mass  and  stiffness 
contributions,  proximity  to  high  stress  areas,  and  the  design  functionality  role.  The 
nodes  of  the  features  found  to  be  non-signihcant  are  removed  from  the  tree  along 
with  the  associated  links.  One  drawback  to  this  approach  is  that  modern  CAD 
systems  rarely  have  the  features  readily  available  by  means  of  a  CSG  tree. 

Date  and  Nishigaki  [9]  present  a  alternative  mesh  simplihcation  method  through 
feature  recognition  and  removal  of  blind  holes,  through  holes,  and  bosses.  The  mesh 
is  segmented  into  regions  enclosed  by  the  edges  of  the  features.  The  edges  are 
extracted  using  the  dihedral  angles  between  common  edge  faces.  Upon  mesh  seg¬ 
mentation,  the  regions  having  an  area  larger  than  the  given  threshold  are  classihed 
as  base  surfaces.  The  triangles  not  within  the  base  surfaces  are  extracted  as  feature 
constructing  triangles  (FCTs).  For  each  FCT  the  feature  type  and  parameters  are 
determined,  the  feature  type  is  identihed  by  the  number  of  boundary  loops  and  its 
concavity,  while  the  dimensions  are  calculated  using  a  least  squares  £t.  Recognized 
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features  are  removed  by  deleting  the  FCTs  and  then  Filing  the  boundary  loops  by 
triangulation.  An  advantage  to  the  proposed  method  is  that  the  removed  features 
are  capable  of  being  recovered  later  on  in  the  process,  which  is  extremely  useful  for 
a  design  model  within  the  industry. 

Gao  et  ah  [14]  proposed  another  framework  for  feature  suppression  based 
CAD  mesh  model  simplihcation.  There  are  three  overall  steps;  segmentation  of  the 
mesh  model,  recognition  of  form  features,  and  feature  suppression. 

Segmentation  of  the  mesh  model  is  achieved  through  an  improved  watershed 
method.  The  first  step  within  the  improved  method  is  to  refine  the  mesh  by  adding 
new  vertices  to  the  model,  which  allow  for  the  model  curvature  to  be  calculated. 
Adding  new  vertices  must  not  change  the  model  geometry,  so  a  different  manner  is 
used  to  add  the  vertices  for  all  the  triangle  edges  as  feature  edges  and  for  triangle 
edges  which  are  not  all  feature  edges.  Upon  rehning  the  mesh,  the  mesh  is  then 
segmented  using  a  multi-descent  strategy.  The  refined  mesh  of  the  segmented  re¬ 
gions  is  then  replaced  with  the  original  mesh  of  the  model,  so  the  original  triangles 
represent  the  segmentation  results.  These  segmentation  results  are  over  segmented 
due  to  the  added  vertices,  therefore  iterative  region  merging  is  used  to  generate  the 
hnal  segmentation. 

After  segmentation  of  the  model  is  complete,  a  graph-based  method  recognizes 
the  features  within  the  model.  Initially  the  region  adjacency  graph  (RAG)  is  set  up, 
where  each  node  refers  to  a  region  and  each  arc  represents  the  relationship  between 
the  regions.  All  base  regions  are  then  determined  from  the  regions  which  have  either 
an  inner  loop  or  concave  vertex  and  deleted  from  the  RAG.  The  remaining  sub- 
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graphs  of  the  RAG  may  then  correspond  to  a  form  feature  capable  of  suppression. 
To  determine  which  ones  are  capable  of  suppression,  the  adjacency  of  the  regions 
with  the  base  features  on  inner  loop  edges  or  concave  edges  are  checked  and  the 
sub-graphs  that  are  adjacent  is  taken  a  s  a  candidate.  Each  feature  graph  candidate 
is  compared  to  the  RAG  of  the  pre-dehned  features  and  a  feature  is  recognized,  if  a 
match  is  found. 

To  generate  the  simplified  model  the  features  recognized  in  the  previous  step 
must  then  deleted  and  the  hole  hlled.  [14]  proposed  three  different  hole  hlling  tech¬ 
niques  planar  Delaunay  triangulation  method  for  planar  features,  poisson  equation 
method  for  curved  features,  and  a  blend  feature  method. 


2.4.2  Blend  Simplification 

Joshi  and  Dutta  [21]  propose  an  algorithm  to  suppress  blends  with  a  circular 
cross-section.  Recognition  of  the  blend  is  determined  from  the  curvature  variation 
across  the  faces,  calculated  along  the  two  directions,  (U,V),  of  the  parameterized 
NURBS  surface,  X(U,V).The  curvatures  at  any  point  are  given  by  Eqns.  2.4  and 
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The  face  is  marked  as  a  potential  blend,  if  Ku  remains  constant  in  the  V 
direction  and  is  greater  than  the  calculated  threshold  and  vice  versa.  The  threshold 
is  calculated  based  upon  the  object  boundary  box,  causing  the  recognition  to  be 
independent  of  the  model  scale.  For  the  potential  blends  to  be  classihed  as  blends, 
the  curves  in  the  direction  of  the  blend  must  be  opposite  an  angle  between  0  and 
180  at  the  spine  curve. 

Blends  are  then  suppressed  in  the  reverse  order  of  creation,  to  maintain  a  valid 
model.  The  blend  order  is  determined  by  visiting  each  blend  and  then  analyzing 
the  adjoining  faces  to  determine  whether  or  not  these  faces  are  also  blends.  If  the 
adjoining  face  is  determined  to  be  a  blend,  an  evaluation  is  done  to  whether  the 
two  blends  were  created  within  the  same  operation  or  if  the  adjoining  faces  blend 
preceded/followed  the  initial  blend  operation.  Suppression  of  the  blend  involves 
removing  the  blend  face  and  then  extrapolating  the  adjacent  faces  to  a  point  of 
intersection,  generating  a  simplihed  model. 

Venkataraman  And  Sohoni  [41]  presented  a  similar  algorithm  for  blend  detec¬ 
tion  and  order  arrangement;  although  rather  than  considering  the  curvatures  over 
the  entire  face  along  the  parametric  direction,  [41]  compares  the  curvatures  of  ad¬ 
joining  faces  at  the  common  edge  midpoints.  Then  the  suppression  algorithm  deletes 
the  blend  faces  and  collapses  the  cross  edges,  mutual  edges  of  the  blend  face  and 
adjacent  edges,  into  vertices  and  the  spring  edges,  the  edges  between  the  blend  face 
and  adjoin  faces,  into  a  single  edge. 

Zhu  and  Menq  et  ah  [42]  also  present  a  method  to  suppress  constant  radius 
rounds  and  hllets,  but  within  a  B-rep  model  focusing  on  the  repairing  the  suppressed 


21 


features.  Initially,  trace  faces,  faces  generated  from  edge  blending,  are  identified  by 
the  convexity  of  their  topological  entities.  These  trace  faces  are  then  removed  from 
the  model,  suppressing  the  rounds  and  hllets.  This  suppression  leaves  gaps  between 
the  adjacent  neighboring  faces  in  the  model.  To  close  this  gap  a  knitting  process  is 
implemented,  which  generates  new  edges  and  vertices.  The  overall  knitting  process 
is  determined  on  the  basis  of  on  two  factors  of  the  H.E  of  the  trace  face,  is  the  H.E. 
a  disk  or  ring  and  is  the  H.E.  a  closed  loop  or  open  loop.  This  process  has  a  couple 
advantages;  the  new  B-rep  model  maintains  geometric  and  topological  consistency, 
along  with  being  a  reversible  process. 

2.4.3  Arbitrary  Shaped  Feature  Simplification 

Joshi  and  Dutta  [21]  proposed  algorithms  for  the  identihcation  and  suppression 
of  different  features  in  B-rep  models.  The  first  algorithm  described  is  for  ’’hole” 
features,  which  are  described  as  any  loop  of  edges  without  a  surface  on  the  inside. 
Initially,  the  edges  having  a  surface  on  one  side  and  no  surface  on  the  other  are 
labeled  as  free  edges.  Upon  labeling  the  free  edges,  all  the  free  edges  are  formed 
into  closed  loops;  starting  at  one  vertex  and  following  the  edge  until  the  starting 
vertex  is  reached  again.  All  these  closed  loops  are  considered  to  be  holes,  even 
the  boundary  of  the  model.  The  boundary  must  be  eliminated  as  a  hole,  so  to 
do  this,  the  perimeters  of  all  the  closed  loops  are  calculated  and  the  hole  with  the 
largest  perimeter  value  is  designated  as  the  model  boundary;  thus,  removing  the 
consideration  for  suppression.  All  loops  with  a  perimeter  value  below  the  given 


22 


threshold  are  deleted  and  replaced  by  a  surface  patch  with  the  same  underlying 
surface  equations  as  the  loop. 

Boss  features,  compound  feature  formed  by  a  planar  surface  and  two  hllet 
rings,  have  also  been  implemented  within  the  framework  designed  by  [21],  For  each 
planar  face  within  the  model,  the  adjacent  faces  on  the  outermost  loop  are  visited 
to  determine  if  they  are  a  blend  (dehned  within  the  blend  section).  If  the  adjacent 
faces  are  found  to  be  blends  and  the  blends  are  all  part  of  the  same  loop,  then 
the  algorithm  analyzes  the  adjacent  faces  of  the  opposing  side  of  the  blends.  If 
the  opposing  faces  are  found  to  either  have  another  blend  loop  on  its  other  side  or 
be  a  blend  loop  itself,  then  the  faces  are  grouped  together  and  labeled  as  a  boss. 
Suppression  of  the  boss  entails  removal  of  all  the  components  and  patching  the 
surface,  similar  to  the  suppression  of  the  holes. 

These  algorithms,  in  combination  with  the  algorithm  proposed  for  blends,  can 
be  extended  to  recognize  and  suppress  any  feature  type  in  a  model. 

2.5  Techniques  based  on  Dimension  Reduction 

This  section  reviews  the  techniques  for  which  the  simplihcation  is  based  upon 
the  dimension  reduction  and  are  organized  into  two  classihcations;  medial  axis  trans¬ 
form  and  mid-surface  abstraction.  Dimensional  reduction  can  be  very  benehcial  for 
FEA.  In  some  cases  it  is  possible  to  reduce  the  dimension  of  the  model  without 
critically  affecting  the  results  of  the  model.  This  dramatically  reduces  the  computa¬ 
tional  time  of  the  model.  For  example,  a  long  round  slender  bar  of  uniform  diameter 


23 


can  be  reduced  to  a  ID  beam  [40].  Other  methods  have  led  to  signihcant  advances 
in  thin  structure  model  simplihcation,  but  may  not  be  applicable  to  generic  struc¬ 
tures;  such  as  an  algebraic  reduction  for  beams  [19],  nonlinear  algebraic  reduction 
for  snap-£t  simulation  [20],  and  an  algebraic  reduction  of  3D  plates  [31],  a  surface 
integration  approach  for  3D  beams  [36]  and  a  method  to  decompose  a  model  into 
parts  and  then  apply  the  mid-surface  abstraction  technique  by  [7]. 

2.5.1  Medial  Axis  Transform 

The  medial  axis  transform,  initially  proposed  by  Blum  [16],  provides  a  skeleton- 
like  representation  of  a  geometric  shape,  which  is  comprised  of  a  set  of  geometric 
entities.  The  geometric  entities  are  generated  by  tracing  out  the  centroid  of  the 
maximal  inscribed  disk/sphere  as  it  rolls  around  the  interior  of  the  object.  This 
concept  provides  both  geometrical  and  topological  proximity  information,  enabling 
a  dimension  reduction  of  the  model,  i.e.  a  2D  model  can  be  reduced  to  a  ID  model 
and  a  3D  model  can  be  reduced  to  a  2D  model. 

Robinson  et  al.  [35]  discusses  a  medial  axis  transform  method  for  both  2D  and 
3D  models.  Only  the  3D  reduction  will  be  touched  on  within,  but  the  paper  can  be 
referenced  for  the  2D  or  mixed  dimension  reduction.  The  first  step  in  the  process 
is  to  identify  the  thin  regions  within  the  model.  In  this  step  medial  flaps,  medial 
edge  or  face  that  terminates  with  a  value  of  zero;  medial  faces  not  defined  by  two 
objects;  and  medial  faces  are  with  a  large  variation  in  3D  medial  radius  or  thickness 
are  all  eliminated  from  the  possibility  of  reduction.  The  remaining  medial  faces  are 
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analyzed  to  determine  if  their  lateral  face  dimensions  are  sufficiently  large  relative 
to  their  thickness.  For  this  a  conservative  aspect  ratio,  since  the  smallest  distance 
is  divided  by  the  largest  thickness,  is  calculated  from  a  2D  MAT  on  the  3D  medial 
face.  A  medial  face  with  an  aspect  ratio  less  than  the  critical  aspect  ratio  is  also 
eliminated  from  the  possibility  of  dimension  reduction.  Not  all  medial  faces  that 
are  left  are  capable  of  reduction,  so  only  the  capable  regions  should  be  reduced. 

After  determining  the  capable  faces  for  reduction,  mid-faces  are  created.  Mid¬ 
faces  are  generated  by  an  insetting  of  the  interface  by  a  distance  equal  to  a  variable 
OS  multiplied  by  the  region  thickness,  achieved  by  reusing  the  2D  MAT  on  each  3D 
medial  face,  and  introducing  rounding  splines  at  convex  corners  of  the  inset.  Upon 
generation  of  the  mid-faces,  the  model  is  partitioned  into  thin  sheet  and  complex 
regions.  The  bounding  edges  of  the  mid-face  is  project  onto  the  medial  faces  and 
the  projected  area  is  deleted  and  replaced  by  a  thin  sheet.  The  thin  sheet  is  then 
extended  to  the  boundary  of  the  component.  Mixed  dimensional  coupling  is  required 
to  at  the  interface  of  the  detailed  section  and  reduced  section  to  join  the  elements. 

The  has  been  extensive  research  done  for  the  MAT  and  only  one  of  the  pro¬ 
cesses  has  been  discussed  within  this  thesis.  Other  algorithms  that  have  been  pre¬ 
sented  are  a  method  for  the  reduction  of  2D  planar  or  3D  shell  models  by  a  Donaghy 
et  al.  [11],  polyhedral  model  homotopy  preserving  medial  axis  transform  by  Sud  et 
al.  [38] ,  and  an  extension  on  the  common  MAT  method  known  as  0-MAT  for  poly¬ 
hedral  mesh  models  by  Foskey  et  al.  [13].  MAT  has  been  extensively  researched  and 
the  survey  paper  by  Attali  et  al.  [5]  can  be  reviewed  for  a  more  complete  review  of 
MAT. 
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2.5.2  Mid-Surface  Abstraction 


Rezayat  [34]  introduced  a  mid-surface  abstraction  technique.  Abstraction  of 
the  mid-surface  involves  four  steps;  pairing  surfaces,  surface  topology  and  size  based 
adjacency  graphs,  mid-surface  patch  generation  by  geometric  interpolation  and  2D 
Boolean  operations,  and  adjacency  based  patch  sewing.  Within  the  surface  pairing 
step,  all  the  faces  of  the  model  are  paired  excluding  the  end-caps  (faces  on  the 
edges  of  the  solid  model)  and  orphans  (faces  that  are  not  on  the  thin- wall  section). 
The  excluded  faces  are  determined  by  whether  or  not  the  ratio  of  the  minimum  face 
height  and  length  divided  by  the  thickness  measure  is  greater  than  a  given  parameter 
X.  A  ratio  less  than  X  distinguishes  a  face  as  either  an  end  face  or  orphan.  To  pair 
the  faces  an  arbitrary  face  is  selected  as  a  seed  face  and  a  ray  is  cast  into  the  material 
side  in  the  normal  direction  of  the  face.  The  face  upon  which  the  ray  hits  is  identihed 
along  with  the  common  edge  faces  as  processed.  This  is  then  repeated  until  all  the 
faces  have  been  processed. 

Upon  pairing  the  faces,  an  adjacency  graph  is  then  generated.  The  nodes  of  the 
graph  are  the  faces  and  the  connecting  arcs  are  the  relationship  between  the  faces, 
either  paired  or  stitched  (common  edge  faces).  The  cyclic  pattern  of  the  adjacency 
graph  references  the  geometric  conhguration  of  the  model.  A  geometric  interpolation 
of  the  paired  faces  is  used  to  generate  a  mid-surface  patch,  where  a  2x2  grid  and 
a  15x15  grid  are  utilized  to  create  the  interpolation  points  of  planar  and  free-form 
surfaces,  respectively.  The  patches  are  sewn  together  at  their  intersecting  edges, 
trimming  the  patch  area  beyond  the  intersection.  If  no  intersection  is  possible,  then 
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a  new  surface  is  created  between  the  patches  by  a  mesh  of  curves  approximation. 
The  part  thickness,  T,  at  a  point,  P,  on  the  mid-surface  of  the  parent  surfaces,  SI 
and  S2,  is  given  by 


T  =  dist  (P,  ^2)  +  dist  (P,  ^1)  (2.6) 

The  authors  have  reported  several  benehts  of  this  method  in  comparison  to 
the  medial  axis  transform  methods.  These  benehts  include  volume  preserving,  ef- 
hcient  feature  suppression  and  detail  removal,  shape  dehnition  through  geometric 
reasoning,  and  part  form  rehection. 

2.6  Sumamry 

The  techniques  described  above  are  present  many  methods  for  which  a  model 
can  be  simplihed  for  different  applications.  For  the  work  discussed  within  this  the¬ 
sis,  the  techniques  based  on  explicit  features  and  dimension  reduction  are  most 
useful;  since  our  main  focus  is  to  suppress  non-critical  features  for  hnite  element 
analysis  purposes.  Within  the  category  of  explicit  features  we  explain  techniques 
for  suppressing  and  modifying  different  features,  such  as:  blind  holes,  through  holes, 
chamfers,  and  rounds  of  any  given  geometric  model.  Although  these  methods  dis¬ 
cuss  ways  in  which  a  model  can  be  simplihed,  the  literature  does  not  provide  error 
estimations  between  the  analyses  of  the  original  and  simplihed  models.  For  this 
[17]  developed  a  formal  theory  for  computing  the  defeaturing-induced  engineering 
analysis  errors,  which  would  help  to  identify  the  suppressible  features.  [23]  presents 
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an  algorithm  or  computing  the  sensitivity  of  the  design  to  hnite  feature  suppression, 
thereby  identifying  the  features  that  would  affect  the  FEA.  The  approach  discussed 
within  this  thesis  has  compared  the  analysis  values,  to  ensure  the  technique  is  within 
an  error  range  of  5%  between  the  models. 


Chapter  3 


Part  Model  Simplification  for  Finite  Element  Analysis 
3.1  Overview 

The  main  focus  of  the  approach  is  to  design  an  automated  tool  that  assists 
a  user  in  identifying  the  features  capable  of  suppression.  The  overall  approach 
towards  generating  a  CAD  simplihcation  tool,  shown  in  Fig.  3.1,  can  be  broken 
down  into  stages.  First,  the  experts’  knowledge  is  captured  through  demonstration 
and  then  suppression  rules  are  generated  with  the  use  of  a  statistical  induction 
learning  technique  [22] .  In  order  to  encode  the  experts’  knowledge  a  representation 
is  needed  with  the  adequate  expressivity.  A  representation  was  developed  that 
allows  the  suppression  rules  to  be  expressed  in  a  neutral  format  and  these  rules  are 
expressed  as  a  binary  decision  tree.  The  developed  language  has  a  very  rich  library, 
which  allows  the  user  to  edit  the  rules  if  necessary.  Upon  learning  the  rules,  they 
are  applied  for  identihcation  of  suppressible  features  within  the  test  parts.  Features 
that  are  found  to  be  suppressible  are  further  examined  based  on  their  parent _child 
relationships.  Lastly,  the  expert  verihes  the  suppressibility  of  the  features  and  the 
CAD  model  is  simplihed  accordingly. 

The  decision  to  suppress  or  preserve  a  feature  strongly  depends  upon  the 
analysis  context  and  the  application.  For  developing  a  rule  language  the  following 
four  attributes  that  play  a  major  role  in  determining  whether  or  not  the  feature  is 
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Verify  Features 
to  be  Suppressed 


Figure  3.1:  Overview  of  the  approach:  Acquire  knowledge  through  expert  demon¬ 


strations,  learn  suppression  rules  from  the  acquired  knowledge  data,  then  apply  the 


learned  rules  for  feature  suppression 


suppressed  were  considered: 

1.  Feature  types  play  a  major  role  in  determining  what  features  should  be  sup¬ 
pressed  and  what  features  should  be  preserved.  Certain  feature  types  are  often 
more  likely  to  be  suppressed.  For  example,  holes,  hllets,  and  chamfers  often 
do  not  signihcantly  affect  the  analysis  results.  However,  the  presence  of  a 
single  small  hole  in  the  model  can  signihcantly  increase  the  size  of  mesh  and 
also  lower  the  quality  of  mesh  element.  So,  such  cosmetic  features  are  prime 
candidates  for  suppression. 
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2.  Feature  dimensions  play  a  major  role  in  determining  whether  they  should  be 
suppressed  or  not.  Often  features  with  small  sizes  are  prime  candidates  for 
suppression.  However,  in  many  cases  it  is  not  the  absolute  dimensions,  but 
the  relative  dimensions  of  the  features  that  determine  whether  they  should  be 
considered  for  suppression  or  not. 

3.  Another  important  consideration  when  selecting  features  for  suppression  is 
their  topological  proximity  to  the  features  on  the  part  on  which  boundary  and 
load  conditions  are  applied.  Features  that  lie  close  to  either  the  point  of  load 
application  or  the  displacement  constraints  should  not  be  eliminated. 

4.  In  CAD  models,  feature  references  can  be  used  to  dehne  placement  of  other 
features.  Such  dependency  among  features  may  cause  problems  during  feature 
suppression.  For  example,  let  us  assume  that  a  feature  is  eligible  for  suppres¬ 
sion.  However,  if  a  large  number  of  features  use  this  feature  as  a  reference, 
then  suppressing  this  feature  may  require  redehning  features  using  other  fea¬ 
tures  as  references.  Not  doing  so  may  lead  to  deletion  of  those  features  as 
well. 

Rules  that  describe  which  features  are  possible  candidates  for  suppression  can 
be  quite  complex  in  nature.  As  more  complex  parts  are  simplihed,  the  rule  library  is 
expanded  upon  in  the  decision  tree.  Conditions  in  the  decision  tree  can  be  expressed 
using  the  values  returned  by  the  following  functions  supported  by  our  representation 
framework: 

Feature.Type  {feature Jd):  This  function  returns  the  type  of  the  feature  identified 
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by  variable  feature Jd. 


Feature_Parameter_Value  {feature-id,  parameter-name):  This  function  returns  the 
value  of  the  feature  parameter  specihed  by  variable  parameter-name.  For  example, 
parameter-name  can  be  set  to  DIAMETER  to  get  the  diameter  of  a  hole  feature. 

Find_Features  (expression):  This  function  returns  the  list  of  features  that  meet 
the  criteria  specihed  by  the  variable  expression.  Expressions  used  as  an  input  for 
this  function  can  be  composed  by  referring  to  features  types,  feature  labels,  and 
keywords.  We  also  support  all  standard  logical  constructs.  For  example,  expression 
((Featnre-Typeithis-feature)  ==  HOLE)  A  (Label  ==  BDUNDARY_CDNDITION) 

A  (DIAMETER  <  10)))returns  all  features  in  the  model  that  are  of  type  HOLE,  do 
not  have  label  as  BOUNDARY  .CONDITION,  and  have  DIAMETER  values  smaller  then  10. 
Values  used  in  the  expression  are  assumed  to  have  the  same  units  as  the  units  being 
used  in  the  CAD  model. 

Feature_Bounding_Box  (feature-id):  This  function  returns  the  bounding  box  of  the 
feature. 

Model_Bounding_Box()  ;This  function  returns  the  bounding  box  of  the  CAD  model. 

Max_Distance_Between_A_Feature_And_Feature_List  (feature-id,  feature-list):  This 
function  returns  the  maximum  distance  of  the  feature  dehned  by  variable  featureTd 
with  the  features  in  the  feature-list. 

Min_Distance_Between_A_Feature_And_Feature_List  (feature-id,  featureJist):  This 
function  returns  the  minimum  distance  of  the  feature  dehned  by  variable  feature-id 
with  the  features  in  the  featureJist. 
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Parent_Features  (featureJd):  This  function  returns  the  list  of  features  that  are 
used  as  a  reference  for  constructing  the  feature  dehned  by  variable  featureJd. 

Child_Features  (featureJd):  This  function  returns  the  list  of  features  that  use  the 
feature  dehned  by  variable  featureJd  as  a  reference. 

The  variable  thisjeature  refers  to  the  feature  being  currently  evaluated.  The 
rule  language  was  generated  based  upon  the  four  attributes  discussed  earlier;  feature 
type,  feature  dimensions,  topological  distance,  and  feature  references.  The  feature 
dimensions  and  topological  distance  are  values  that  are  relative  to  the  overall  size 
of  the  3D  model.  To  acquire  a  representative  rule  base  for  3D  models,  these  values 
must  be  normalized  to  eliminate  the  correlation  between  the  overall  model  size  and 
the  feature  values.  To  accomplish  this  task  the  feature  values  were  divided  by  the 
average  of  the  sides  of  the  overall  models  bounding  box.  This  normalizes  the  data 
and  removes  the  dimensional  bias  of  the  of  the  feature  value  with  respect  to  the 
overall  model  size. 

The  representation  allows  a  wide  variety  of  knowledge  associated  with  feature 
suppression  in  a  wide  variety  of  contexts  to  be  easily  captured.  A  decision  tree  in 
the  context  of  brackets  and  mounts  to  illustrate  how  the  approach  works  has  been 
developed.  Section  3.2  describes  the  framework  to  learn  the  rules  for  suppression 
based  on  the  knowledge/data  collected  during  the  expert  demonstrations.  Section 
3.3  presents  an  algorithm  for  applying  these  learned  rules  and  Section  3.4  presents 
test  results  of  the  simplihcation  process. 
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3.2  Using  Inductive  Decision  Tree  to  Generate  Suppression  Rules 


Machine  learning  mechanizes  the  acquisition  of  knowledge  from  experience 
through  the  use  of  computational  methods  [24],  It  aims  to  improve  performance 
by  automating  manual  time-consuming  activity  within  the  knowledge  process  and 
increase  accuracy  by  discovering  and  exploiting  regularities  in  training  data  [24] .  A 
popular  technique  in  robotics  is  learning  from  demonstration  [4] .  Inductive  decision 
trees  have  proven  to  be  a  successful  tool  within  the  numerous  fields  and  these  are 
proven  to  be  advantageous  in  improving  the  accuracy  or  efficiency  of  the  system. 
Many  real  world  problems  have  implemented  the  technique  and  their  experiments 
have  shown  that  the  learned  rules  are  more  accurate  than  the  handcrafted  ones 
previously  used  [15,  12], 

The  framework  implemented  here  requires  the  rules  to  classify  the  features  as 
suppressible  or  non-suppressible.  The  feature  attributes  used  for  the  rule  learning 
process  are  both  discrete,  having  a  predehned  set  of  options,  or  continuous,  real 
numbers.  Regression  and  other  model-htting  methods  require  discrete  attributes  to 
be  encoded  as  numbers,  in  which  the  transformation  introduce  spatial  connotations 
that  may  affect  the  outcome.  On  the  other  hand  Bayesian  methods  require  all 
attributes  to  be  discrete,  so  continuous  attributes  must  be  divided  into  sub-ranges 
[32],  For  these  reasons  decision  trees  were  used,  since  it  can  cope  with  both  discrete 
and  continuous  attributes.  Trees  are  also  very  clear  and  concise  in  style;  along  as 
context  sensitive,  allowing  for  the  relevance  of  different  attributes  to  be  conditional 
on  the  outcomes  of  previous  tests  [32], 
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Figure  3.2:  A  Simple  decision  tree:  An  example  of  a  representative  decision  tree 
for  suppressing  a  hole  feature  considering  attributes  diameter,  minimum  distance 
to  the  boundary  condition  feature,  and  the  geometric  type.  Rectangular  boxes  and 
elliptical  boxes  represent  the  non-terminal  nodes  and  the  leaf  nodes  respectively. 

The  decision  tree  is  used  to  classify  the  suppressibility  of  the  features  in  the 
following  manner.  Given  a  feature  to  classify,  the  decision  tree  is  traversed  from 
the  root  node  to  a  leaf  node  along  a  path.  All  non-terminal  nodes  in  the  decision 
tree  only  define  test  conditions.  These  nodes  in  the  decision  tree  point  to  exactly 
two  children  nodes.  The  left  node  is  the  node  that  should  be  executed  if  the  con¬ 
dition  is  true.  The  right  node  is  the  node  that  should  be  executed  if  the  condition 
is  false.  Terminal  nodes,  leaf  node,  in  the  decision  tree  assign  the  variable  sup- 
pressibihty_state  for  the  feature.  If  the  feature  is  suppressible,  then  the  variable 
suppressibility ^tate  is  set  to  TRUE.  Otherwise,  this  label  is  set  to  FALSE.  For  a 
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given  set  of  rules  as  shown  below,  a  simple  decision  tree  can  be  formed  (Fig.  3.2) 
by  using  various  split  conditions.  The  decision  to  suppress  or  not  to  suppress  can 
be  made  based  on  the  attribute  values.  Features  that  are  found  to  be  suitable  for 
suppression  are  further  examined  based  on  their  parent-child  relationship. 

•  If ( (Feature.Type ( this_feature) ==H0LE)  A  (Feature_Parameter_Value 
( th'is_/eati4re, DIAMETER)  <  0.03)  A 

( (Min_Distnace_Between_A_Feature_And_Feature_List 
(this_feature, FindJFeatnres  (LabeI==BOUNDARY .CONDITION) ) ) 

>=  0.05)  then  (SuppressibiIity_State  =  TRUE)) 

•  If ( (Feature.Type ( this^feature) ==H0LE)  A  (Feature_Parameter_VaIue 
( th'is./eatiire, DIAMETER) >=0 . 03)  A  (Feature_Parameter_VaIue 
(this_f eature, TYPE) ^  BLIND)  A (Feature_Paraineter_Value 

( thisjeature, TYPE)  ^  THROUGH)  then 
(SuppressibiIity_State  =  TRUE)) 

The  set  of  attributes  used  for  the  learning  process  will  be  dependent  upon  the 
feature  type.  Holes  have  the  feature  parameters  of  diameter,  type  (through/blind), 
and  minimum  distance  to  the  boundary  conditions.  The  diameter  and  minimum 
distance  to  the  boundary  condition  attributes  are  both  continuous  in  nature,  whereas 
the  type  attribute  is  discrete  in  nature.  Rounds  have  the  attribute  types  of  radius, 
conic  value  and  minimum  distance  to  the  boundary  conditions.  Chamfers  have  the 
attribute  types  of  Dl,  D2,  and  minimum  distance  to  the  boundary  conditions.  Even 
though  the  attributes  for  rounds  and  chamfers  are  all  considered  to  be  continuous 
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in  nature,  a  decision  tree  is  still  used  for  analysis  consistency  with  the  hole  features. 

These  attributes  were  selected  to  convey  the  key  aspects  of  the  features;  feature 
type,  feature  size,  and  minimum  geometric  distance  to  the  boundary  conditions. 
These  aspects  have  a  effect  on  whether  the  feature  is  critical  or  not.  A  couple 
models  have  been  analyzed  to  demonstrate  the  effects  of  the  different  aspects  on 
the  analysis  results.  The  importance  of  knowing  the  minimum  distance  to  the 
boundary  conditions  is  illustrated  in  Figs.3.3(a)and(b).  There  is  a  round,  meeting 
the  dimensional  requirements  for  suppression,  but  it  is  within  a  critical  distance 
to  the  boundary  condition.  Suppression  of  this  round  feature  is  detrimental  to  the 
accuracy  of  the  results.  Upon  suppressing  this  round  feature  the  critical  stress  value 
changes  from  a  value  of  42.27  MPa  within  the  original  model  illustrated  in  3.3(a)  to 
a  value  of  32.87  MPa  shown  in  3.3(b).  When  the  round  is  maintained  the  critical 
stress  value  is  comparable  to  the  original  model  at  a  value  of  43.61  MPa.  Not  only 
is  value  of  the  critical  stress  erroneous,  but  the  location  of  the  critical  stress  value 
has  changed  and  can  also  be  seen  in  the  following  hgures. 

To  learn  the  decision  trees  the  analysts  experiential  knowledge  was  obtained 
by  conducting  numerous  simplihcation  demonstrations  on  various  training  parts, 
where  the  information  about  the  features  and  their  associated  suppressibility  states 
assigned  by  the  analyst  was  collected.  Collecting  the  data  of  the  non-suppressed 
features  along  with  the  suppressible  features  allows  for  a  more  expansive  rule  base 
for  easy  classihcation. 
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(b)  Simplified  model  witfi  critical  round  suppressed 


Figure  3.3:  Example  of  suppression  based  on  inter-feature  distances 


3.2.1  Training  Parts 

Formation  of  reasonably  good  decision  trees  can  be  possible  by  generating  a 
large  number  of  training  parts  for  demonstration  purposes;  although  this  can  be 
manually  time  consuming  for  the  expert  to  conduct  demonstrations  for  such  a  large 
data  set.  However,  with  the  most  informative  training  examples  the  decision  tree 
can  classify  the  features  correctly  without  a  need  of  as  many  examples,  through 
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Figure  3.4:  Training  examples  used  to  collect  attributes  of  holes,  rounds,  and  cham¬ 
fers.  Synthetic  parts  are  generated  from  these  examples  to  tune  the  transition 
boundaries  of  the  decision  trees  for  better  classification. 

the  use  of  synthetic  parts  to  refine  the  transition  point  values  at  the  classification 
boundaries.  For  this,  some  heuristics  were  used  for  choosing  places  where  the  deci¬ 
sion  tree  performs  poorly  and  generate  an  additional  set  of  training  data  at  those 
places  for  refining  the  boundaries. 

Initially,  10  training  parts  were  started  with  that  offer  diversity  in  terms  of 
feature  types  and  feature  parameters.  A  decision  tree  is  formed  from  these  initial 
parts.  This  decision  tree  failed  to  give  the  correct  answer  for  those  places  where  the 
transition  region  between  the  suppressible  and  non-suppressible  class  is  very  large 
and  hence  locating  the  correct  transition  boundary  separating  these  two  classes  is  a 
difficult  task.  To  overcome  this  issue  the  divide-and-conquer  strategy  (binary  search) 
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Feature  Type 

No.  of 

of  Features 

No.  of 

Suppressible  Seatures 

No.  of 

Non-Suppressible  Features 

Holes 

169 

60 

109 

Chamfers 

57 

23 

34 

Rounds 

118 

75 

43 

Table  3.1:  Training  data  collected  from  the  examples  shown  in  Fig.  3.4. 


was  used  to  locate  this  region  and  then  tune  the  decision  tree  with  the  synthetic 
test  cases  generated  in  the  vicinity  of  the  boundary.  To  determine  the  transition 
region  another  37  synthetic  parts  were  generated  to  rehne  the  transition  points  of 
the  rule  base.  From  these  37  synthetic  parts  the  only  recorded  features  were  those 
features  used  for  the  transition  region  or  features  that  required  alteration  when 
designing  the  transition  feature  of  the  synthetic  part.  Only  recording  these  features 
within  the  synthetic  parts  eliminated  the  chance  of  biasing  the  data,  through  the 
introduction  of  the  same  features  numerous  times  within  the  data.  Since  the  data 
is  being  recorded  in  this  fashion,  there  will  be  a  relatively  low  total  feature  count 
for  the  number  of  total  training  parts.  The  training  parts  can  be  seen  in  Figs.  3.4. 
Numerous  hole,  round,  and  chamfer  features  were  obtained  from  these  examples  for 
rule  learning.  There  was  a  total  of  344  features  within  the  training  parts,  158  were 
suppressible  and  186  were  non-suppressible.  Table  3.1  displays  the  statistics  of  these 
features  within  the  training  parts. 
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3.2.2  Algorithm  for  Generating  Rules 


In  this  work,  a  decision  trees  was  bnilt  from  the  geometric  feature  attributes 
and  the  associated  suppressibility  values.  From  the  set  of  training  parts  a  decision 
tree  [22]  can  be  formed  to  generate  a  classihcation  model  that  can  generalize  beyond 
the  training  samples  so  that  new  samples  can  also  be  classified  with  a  high  accuracy. 
In  order  to  do  this,  the  decision  tree  must  capture  some  meaningful  relationships 
between  the  output  class  (C  =  Gi,  G2, ...,  G^)  and  its  input  data  space  consists  of 
values  from  independent  parameters  (h  =  hi,  ^2,  ...h^).  Within  the  data  there  are 
three  independent  parameters;  for  instance  a  hole  feature  utilizes  the  diameter,  type, 
and  minimum  distance  to  a  boundary  condition  feature.  The  dependent  parameter, 
Gj,  represents  the  class,  whether  suppressible  or  non-suppressible. 

Development  of  the  tree  from  the  training  cases  begins  at  a  test  node,  known 
as  the  root  node,  of  the  tree.  The  tree  builds  from  the  root  node,  splitting  the 
dataset  of  one  particular  attribute  at  each  sub-node.  The  division  of  the  attributes 
is  determined  based  on  the  information  gain,  which  is  maximized  at  each  node;  the 
attribute  with  the  highest  information  gain  is  placed  at  the  node.  The  information 
gain  is  the  entropy  relationship  of  the  node  before  and  after  the  split  based  on  the 
specihed  attribute  [22].  This  process  is  continued  until  the  node  contains  only  one 
class;  then  the  node  is  no  longer  a  test  node,  but  considered  to  be  a  leaf  node.  The 
tree  is  now  considered  to  be  fully  grown.  Fully  grown  trees  tend  to  over-£t  the  data, 
meaning  that  the  tree  characterizes  too  much  detail  or  noise  from  within  the  training 
set.  To  reduce  the  error  of  over-fitting  the  data,  the  tree  was  pruned.  Pruning  the 
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tree  involves  calculating  the  error  rates  of  the  branches  of  the  tree.  Comparison  of 
the  error  rates  determines  whether  the  branch  is  left  unchanged,  changed  to  a  leaf 
node,  or  replaced  by  an  alternate  sub-tree. 

The  generated  decision  trees  for  the  different  feature  types  are  then  written  in 
the  framework  rule  language.  Tables  3. 2-3. 5  display  the  actual  rules  produced  from 
the  decision  trees  of  the  hole,  round,  and  chamfer  features. 

Rule  Suppressibility.State 


•  If  (  (Feature.Type  ( this.feature)  ==  HOLE)  A  (Min_Distnace3etween_A_Feature_And_Feature_List 
(thts_/eatiire,Find_Features  (Label  ==  BDUNDARY.CONDITIQN) ) )  <=  0.02368) 


FALSE 


•  If  (  (Feature.Type  ( this.feature)  ==  HOLE)  A  (Min_Distnace3etween_A_Feature_And_Feature_List 
(thts_/eatiire,Find_Features(Label  ==  BDUUDARY-CONDITION) ) )  >  0.02368) 

A  (Feature_Parameter_Value  (ihis_/eaiure, DIAMETER)  <=  0.08451)  A 
(Feature_Parameter_Value  ( i?iis_/eaiure, TYPE)  ==  THROUGH) 


TRUE 


•  If  (  (Feature_Type  ( this^feature)  ==  HOLE)  A  (Min_Distnace3etween_A3eature_And3eature_List 
(t?iis_/eatiire,Find_Features  (Label  ==  BDUNDARY_C0NDITIDN) ) )  >  0.02368) 

A  (Feature_Parameter_Value(tfiis_featiire,DIAMETER)  <=  0.08451)  A 

TRUE 

(Feature_Parameter_Value  ( this.feature ,  TYPE)  ==  BLIND)  A  (Feature3aramet er.Value 
(this./eature, DIAMETER)  <=  0.018626) 


Table  3.2:  Rules  for  computing  the  suppressibility_state  of  holes 
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Rule 


Suppressibility.State 


•  If  (  (Feature.Type  ( this.feature)--  HOLE)  A  (Min_Distnace_Between_A_Feature_And_Feature_List 
( this.feature , Find_Features  (Label==BOUNDARY_CQNDITIDN)  )  )  >  0 . 02368) 

A  (Feature_Paraineter .Value  (i/iis./eaiure, DIAMETER)  <—  0.08451)  A  FALSE 

(Feature_Parameter_Value  (iAis_/eaiure, TYPE)==BLIND)  A  (Feature Jarameter-Value 
(t?its_/eatiire, DIAMETER)  >  0.018626)  A  (Feature J’arameter.Value 
(t?iis_/eatiire, DIAMETER)  <=  0.02314) 


•  If  (  (Feature.Type  ( this^feature)  ==  HOLE)  A  (Min_Distnace.BetweenJlJ^eature JlndJ^eature_List 
(t;its_/eatiire,Find_Features(Label  ==  BDUNDARY.CONDITION) ) )  >  0.02368) 

A  (Feature_Paraineter .Value  (iAis./eaiure, DIAMETER)  <—  0.08451)  A  TRUE 

(Feature_Parameter .Value  (iAis./eaiure, TYPE)  ==  BLIND)  A  (Feature_Parameter .Value 
(t/iis./eature, DIAMETER)  >  0.018626)  A  (Feature Jarameter.Value 
(t/iis./eature, DIAMETER)  >  0.02314) 


•  If  (  (Feature.Type  ( this.feature')  ==  HOLE)  A  (Min_Distnace3etween_ATeature_And_Feature_List 
(t;its_/eatiire,Find_Features(Label  ==  BOUNDARY-CONDITION)))  >  0.02368) 

FALSE 

A  (Feature-Parameter .Value  (iAis./eaiure, DIAMETER)  >  0.08451)  A 
(Feature-Parameter .Value  ( this. feature, TYPE)  ==  THROUGH) 


•  If  (  (Feature.Type  ( this.feature)  ==  HOLE)  A  (Min_Distnace3etweenA3eatureAnd3eature_List 
(t;its_/eatiire,Find_Features(Label  ==  BOUNDARY-CONDITION)))  >  0.02368) 

A  (Feature-Parameter .Value  (iAis./eaiure, DIAMETER)  >  0.08451)  A 
(Feature-Parameter .Value  ( this. feature, TYPE)  ==  BLIND) 

TRUE 

Table  3.3:  Rules  for  computing  the  suppressibility .state  of  holes 


43 


Rule 


Suppressibility  .State 


•  If  (  (Feature-Type  ( t/its_/eatiire)  ==  CHAMFER)  A 

FALSE 

(Min_Distnace_Between_A_Feature_And_Feature_List  ( fhis_/eaiure,Find_Features(Label 
==  BOUNDARY-CONDITION)))  <=  0.004121) 


•  If  (  (Feature-Type  ( t?i-is_/eatiire)  ==  CHAMFER)  A 

(Min_Distnace3etween_A_Feature_And_Feature_List  ( tftis_/eaiure, Find-Features  (Label 

TRUE 

==  BOUNDARY-CONDITION)))  >  0.004121)  A  (Feature-Parameter-Value  (.this.feature,Dl') 

<=  0.01714) 


•  If  (  (Feature-Type  ( t/its-/eature)  ==  CHAMFER)  A 

(Min-Distnace3etween-A-Feature-And-Feature-List  (  fhis-/eaiure, Find-Features  (Label 
==  BOUNDARY-CONDITION)))  >  0.004121)  A  (Feature3arameter-Value  (.this.feature,Dl')  > 
0.01714) 


FALSE 


Table  3.4:  Rules  for  computing  the  suppressibility .state  of  chamfers 
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Rule 


Suppressibility  .State 


•  If  (  (Feature.Type  ( t?i-is_/eatiire)  ==  ROUND)  A  (Feature_Parameter_Value 

FALSE 

(t?its_/eatiire, RADIUS)  <=  0.02623)  A  (Min_Distnace.Between_A J^eature_And_Feature_List 
(t?its_/eatiire,Find_Features(Label==BOUNDARY_CDNDITIDN)))  <=  0.013545) 


•  If  (  (Feature.Type  ( t/its_/eatiire)  ==  ROUND)  A  (Feature_Parameter_Value 

TRUE 

(t/its_/eatiire, RADIUS)  <=  0.02623)  A  (Min_Distnace3etween_A_Feature_And_Feature_List 
( this.feature , Find_Features  (Label==BOUNDARY_CONDITIDN)  )  )  >  0 .013545) 


•  ((Feature.TypeCtftis./eature)  ==  ROUND)  A  (Feature_Parameter_Value 

FALSE 

(t?its_/eatiire, RADIUS)  >  0.02623)) 


Table  3.5:  Rules  for  computing  the  suppressibility .state  of  rounds 
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3.3  Applying  Learned-Rules  for  Simplification 


The  rules  have  been  learned  from  training  parts,  now  the  next  step  within 
the  approach,  as  seen  in  Fig.  3.1,  is  to  apply  the  rules  to  new  models.  To  do  this, 
a  utility  was  needed  to  extract  the  attribute  information  and  classify  the  features 
within  the  model,  for  this  ProE  was  used. 

3.3.1  Extracting  CAD  Model  Information 

ProE  has  a  Pro/Toolkit  API,  which  allows  customization  within  the  frame¬ 
work.  In  this  API,  the  design  model  is  stored  by  ProMdl,  a  data  handle  to  retrieve 
the  complete  information  of  the  model  such  as  features,  their  dependencies,  etc.  By 
utilizing  the  ProMdl  the  framework  traverses  through  all  the  features,  where  each 
feature  can  contain  a  large  amount  of  information.  In  ProE,  the  complete  design 
model  is  stored  as  a  feature  element  tree  that  contains  all  information  about  the 
feature,  including,  for  example:  all  attributes,  such  as  the  material,  placements  of 
holes;  all  references  to  the  existing  geometries  such  as  placement  references;  and  all 
dimension  values.  For  example,  a  “hole”  feature  contains  information  as  to  whether 
the  hole  is  straight,  standard,  or  sketched;  whether  the  hole  is  a  clearance  hole  or 
a  threaded  hole;  parameters  of  the  hole,  etc;  all  this  information  can  be  obtained 
through  the  element  tree. 

To  implement  the  part  of  the  rule  Feature_Type(/eaf'wre_id),  the  element  tree  of 
the  complete  design  model  was  traversed  and  determined  the  feature  /  corresponding 
to  the  featuredd.  Now,  the  element  tree  of  this  feature  /  was  traversed  to  extract 
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(a)  (D  X  D)  (b)  (D1  x  D2)  (c)  (ANG  x  D) 

.15  H  — 


(d)  (45  X  D)  (e)  (O  X  O)  (f)  (01  x  02) 

Figure  3.5:  Suppression  based  on  chamfer  geometric  types:  (a)  D  x  D  (b)  D1  x  D2 
(c)  ANG  X  D  (d)  45  X  D  (e)  O  X  O  and  (f)  01  x  02 

the  internal  structure  of  the  feature  information,  such  as  the  diameter  value  for  a 
hole  or  the  D1  value  for  a  chamfer.  There  are  many  types  of  features  that  can  be 
present  within  a  GAD  model,  including  holes,  rounds,  chamfers,  hllets,  steps,  and 
slots.  Gurrently,  the  framework  implemented  considers  (a)  chamfers  (b)  holes  and 
(c)  rounds.  In  the  case  of  chamfers  there  is  an  option  of  six  different  geometric 
representations;  described  by  one  distance  (D  x  D)  (Fig.  3.5) (a);  by  two  distances 
(D1  X  D2)  (Fig.  3.5)(b);  by  an  angle  and  a  distance  (ANG  x  D)  (Fig.  3.5)(c);  by 
a  45°  angle  and  a  distance  (45  x  D)  (Figure  3.5) (d);  by  an  off-set  distance  (O  x 
O)  (Fig.  3.5) (e);  or  by  two  off-set  distances  (01  x  02)  (Fig.  3.5)(f).  In  a  similar 
manner  the  suppressibility^tate  of  holes  is  described  by  two  types;  blind  holes  and 
through  holes;  and  rounds  are  described  as  either  conical  or  circular. 

In  addition  to  the  size  and  type,  feature  locations  must  also  be  calculated 
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for  the  rule  implementation.  Feature  locations  play  an  important  role  in  determin¬ 
ing  the  extent  of  the  feature  criticality,  since  often  features  that  are  very  close  to 
the  boundary  conditions  are  very  critical  and  should  not  be  eliminated  irrespec¬ 
tive  of  their  size.  In  most  situations,  closeness  in  terms  of  the  topological  distance 
is  more  important  than  the  geometric  distances.  Features  that  are  far  away  from 
the  boundary  conditions  can  also  become  critical  if  they  are  located  within  the 
stress  flow  path  from  the  load  application  to  the  displacement  constraints.  Utilizing 
the  available  ProE  API  functions  each  and  every  feature  of  the  design  model  was 
visited,  computing  the  inter-feature  distances  among  all  the  relevant  feature  pairs 
available  in  the  model.  Then  the  feature  pairs  are  retrieved  associated  with  load  or 
displacement  constraints  locations.  In  the  framework,  features  associated  with  load 
or  displacement  constraint  locations  are  labeled  as  BOUNDARY  .CONDITION  features. 
This  labeling  needs  to  be  done  manually. 

3.3.2  Utilizing  CAD  Model  Information  within  Rules 

The  collected  attribute  data  for  each  feature  in  the  model  from  section  3.3.1 
was  normalized  just  like  the  training  data,  to  remove  the  dimensional  bias  of  the 
parts.  The  normalized  data  was  then  evaluated  utilizing  the  machine  learned  rules, 
generated  from  the  training  data.  For  instance,  the  attribute  values  of  a  hole,  con¬ 
sisting  of  the  diameter,  type,  and  minimum  distance  to  the  boundary  condition, 
collected  in  section  3.3.1  were  compared  to  the  rule  values  in  Table  3.2.  Referring 
to  Table  3.2  a  hole  with  a  calculated  normalized  minimum  distance  to  the  boundary 
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condition  features  of  0  would  be  classified  as  non-suppressible  without  any  other 
attribute  analysis,  where  as  a  normalized  distance  of  .03  would  then  cause  the  di¬ 
ameter  attribute  to  be  analyzed.  Upon  analyzing  the  diameter,  the  type  would  then 
be  used  for  comparison  and  depending  on  the  diameter  and  type  values  the  analysis 
could  potentially  stop  there  or  continue.  After  the  feature’s  attribute  values  were 
analyzed  using  the  generated  rule  set,  the  features  were  then  labeled  as  suppressible 
or  non-suppressible  determined  by  the  classification  obtained  from  the  rule  base. 


Actual 

Suppressibility. State 

Rule 

Suppressibility  .State 

Classification 

Error  Type 

True 

True 

Correct 

None 

False 

False 

Correct 

None 

True 

False 

Incorrect 

Conservative 

False 

True 

Incorrect 

N  on-  Conservative 

Table  3.6:  Errors  and  error  types  within  the  rule  learning  applieation:  Conserva¬ 
tive  errors  can  be  allowed  as  they  will  not  affect  the  analysis  results  but  Non- 
Conservative  errors  can  generate  erroneous  results. 


Implementation  of  the  learned  rules  from  the  training  data  may  sometimes 
be  flawed,  which  can  generate  two  different  error  types,  shown  in  Table  3.6  when 
compared  with  the  actual  suppressibility .states  of  the  features.  The  hrst  type  is  a 
conservative  error,  which  is  an  error  that  will  not  affect  the  analysis  results  of  the 
simplihed  model.  This  occurs  when  the  rule  states  not  to  suppress  a  feature  that  has 
an  actual  suppressibility .state  of  TRUE,  or  non-critical  feature.  The  other  type  is  a 
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non- conservative  error,  which  will  likely  affect  the  analysis  results  of  the  simplified 
model.  This  error  occurs  when  the  rule  states  to  suppress  a  feature  that  has  an 
actual  suppressibility -Status  of  FALSE,  or  critical  feature.  For  the  reason  that  non¬ 
conservative  errors  can  generate  erroneous  results,  the  error  has  been  limited  within 
the  inductive  learning  process  to  allow  for  only  conservatives  errors.  This  was  done 
through  the  use  of  the  synthetic  parts  to  refine  the  transition  points. 

3.3.3  Accounting  for  Parent-Child  Relationships  during  Feature  Sup¬ 
pression 

Design  model  generation  using  a  feature-based  approach  allows  flexibility  to 
compile  a  set  of  information  about  the  product  at  higher  levels  of  abstractions.  This 
information  can  be  represented  in  the  form  of  data  structures  that  maintain  the  asso¬ 
ciations  between  the  geometric  entities,  further  aiding  in  making  any  changes  to  the 
design.  All  major  parametric  CAD  modeling  systems  have  adopted  the  feature-based 
design  approach  by  representing  the  final  design  model  with  features  dependent  on 
one  another.  This  is  called  the  design  feature  history  tree  [33,  37].  The  same  design 
model  may  be  designed  in  several  different  ways,  depending  on  the  designer,  who 
builds  the  design  feature  history.  For  the  finite  element  analysis  the  order  in  which 
one  simplifies  the  features  is  immaterial.  However,  it  is  necessary  to  understand 
the  dependent  relationships  among  the  features,  as  the  suppression  of  one  feature 
may  suppress  another  feature  that  may  be  very  critical  to  the  analysis.  Hence,  it  is 
imperative  to  understand  features’  parent-child  relationships  for  suppressions.  For 


50 


example,  on  a  block  a  hole  was  first  created  and  then  a  round  operation  was  per¬ 
formed  on  it.  The  round  is  dehned  as  the  child  of  the  parent  hole,  which  in  turn  is 
the  child  of  the  parent  block.  In  this  case,  suppressing  the  hole  would  also  suppress 
the  round.  Extracting  these  parent-child  dependencies  is  crucial  in  order  to  avoid 
suppressing  any  non-critical  features. 

In  the  approach,  the  feature  dependencies  were  initially  extracted  from  the 
model  using  the  available  API  functions  and  represented  as  a  feature  dependency 
graph,  as  illustrated  in  Figure  3.6.  Let  G{V,  E)  represent  the  parent-child  relation¬ 
ship  graph.  Each  vinV  represents  a  feature  in  the  feature  tree.  Each  element  e 
in  E  is  a  directed  edge  (n,  v)  that  represents  the  fact  that  node  n  is  a  parent  of 
node  V.  Since  features  are  inserted  one  by  one  in  the  feature  tree  and  can  only 
use  the  already  dehned  features  as  references,  by  the  property  of  the  feature  tree 
construction,  G  is  a  directed  acyclic  graph. 

Within  the  graph  the  nodes  represent  the  various  features  of  the  model.  Each 
dependency  is  represented  by  an  arrow,  which  is  directed  from  the  parent  feature 
to  the  child  feature.  The  tree  illustrates  the  dependencies  for  easy  analysis.  Af¬ 
ter  obtaining  the  G  value,  we  follow  an  algorithm  for  identifying  different  cases  of 
suppression  possibilities,  as  a  feature  should  not  be  suppressed  just  because  it  is 
suppressible.  It  may  have  child  features  indicating  that  it  is  not  suppressible.  On 
the  other  hand,  it  can  be  suppressed  if  all  of  its  direct  children  can  be  suppressed. 
Hence,  we  needed  to  use  a  special  sequence  to  determine  which  features  can  be 
suppressed  by  accounting  for  their  parent-child  relationships. 

The  algorithm  for  suppressing  features  is  as  follow: 
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1.  Construct  G  by  inserting  a  node  in  V  corresponding  to  every  feature  in  the 
feature  tree.  If  feature  /'  uses  a  feature  /  as  a  reference,  then  insert  an  edge 
(n,  /')  in  G.  Here  vertex  v  corresponds  to  feature  /  and  vertex  v'  corresponds 
to  feature  /'. 

2.  Find  nodes  X  in  V  that  do  not  have  parents  and  have  label  suppressibil- 
ity_state  =  TRUE.  If  X  is  empty,  then  reset  suppressibihty_state  for  all  nodes 
to  be  FALSE  and  stop.  Otherwise  do  the  following: 

(a)  Remove  every  node  in  X  from  G  by  deleting  associated  vertex  and  edges. 

(b)  Go  to  Step  2. 

Figure  3.6(a)  shows  a  parent-child  relationship  graph.  Figures  3.6(b)  through  (e) 
show  how  executing  the  above  algorithm  correctly  identihes  the  features  to  be  sup¬ 
pressed.  This  algorithm  discovers  at  least  one  feature  to  suppress;  otherwise  it 
terminates  in  Step  2.  Hence,  it  takes  at  most  V  iterations  to  complete  its  task.  In 
this  way,  what  our  tool  is  capable  of  suppressing  the  features  based  on  the  feasibility 
dependencies. 
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Features  with  suppressability  status 
set  to  TRUE 


(a)  Initial  tree  (b)  First  pass  for  suppression 


(c)  Second  pass  for  suppression 


(d)  Third  pass  for  suppression 


(e)  Fourth  pass  for  suppression 

Figure  3.6:  Dependency  graph  showing  an  example  of  finding  features  to  be  sup¬ 
pressed  based  on  the  parent-child  relationships.  Dependent  relationships  among  the 
features  should  be  considered  before  suppression;  the  suppression  of  one  feature  may 
suppress  the  dependent  critical  feature. 


53 


3.4  Results  &  Discussion 


To  demonstrate  the  simplification  tool  five  parts  were  considered  Fig.3.7- 
Fig.3.11.  The  simplified  models  were  obtained  by  using  the  learned  rules  from  the 
training  models.  Now  a  comparison  was  done  of  the  analysis  results  of  the  original 
and  simplified  models  (Fig.  3.7(b)  -  Fig. 3. 11(b).  In  terms  of  analysis  we  initially 
add  displacement  constraints  and  applied  forces  were  added  to  each  part  and  a  fi¬ 
nite  element  analysis  was  conducted  with  the  help  of  Pro/Mechanica  (available  in 
ProE).  Pro/Mechanica  requires  every  model  to  contain  a  constraint  and  a  force  for 
analysis. 

The  location  of  the  displacement  constraints  along  with  the  applied  forces 
is  explained  below.  Figs.  3.7(a)  -  Fig.3. 11(a)  illustrate  the  boundary  conditions 
(constraints  and  forces)  applied  to  the  to  the  five  complex  models. 

1.  Part- A;  The  constraints  were  placed  at  the  eight  outside  holes  on  the  base 
of  the  part,  where  it  is  bolted  into  place.  The  forces  were  applied  in  the 
downward  vertical  direction  located  at  the  two  center  holes  of  the  extension 
at  the  top  of  the  part  (Fig.  3.7(a)). 

2.  Part-B:  The  constraints  are  located  at  the  vertical  slots.  A  downward  force 
was  applied  on  the  pin  holes  shown  in  Fig.  3.8(a). 

3.  Part-C:  The  constraints  are  located  at  the  rectangular  slots  in  the  model  that 
are  connection  points  to  a  vertical  piece.  A  distributed  vertical  load  was  placed 
on  the  plane  of  the  face  perpendicular  to  the  vertical  plane  (Fig.  3.9(a)). 
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4.  Part-D;  The  constraint  was  applied  at  the  location  of  the  hole,  where  this  part 
will  be  bolted  into  place.  Mnltiple  forces  were  applied  at  different  locations 
of  the  model,  (1)  tensile  forces  were  applied  at  the  hole  on  the  free  end  and 
(2)  another  force  was  applied  on  the  other  end  of  the  part  in  the  opposite 
direction  at  the  middle  pin  shown  in  Fig.  3.10(a). 

5.  Part-E;  The  part  slides  into  grooves  on  a  wall  snpport  and  then  pinned  into 
place.  The  top  cnrve  latches  on  a  lip  and  acts  as  the  constraint  of  the  model.  A 
force  was  applied  in  the  vertical  direction  on  the  horizontal  face  (Fig.  3.11(a)) 
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(a)  Boundary  conditions 


Hole  features 
(4  comers) 
suppressed 


Round 
features 
(witliiii  slots ) 
suppressed 


(b)  Suppression  of  features 


Figure  3.7:  Part  A:  boundary  condition  features  and  feature  suppression  after  rule 


implement  at  ion 
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(a)  Boundary  conditions 


Chamfer  features 
arc  suppressed 


Chamfer  features 
arc  suppressed 


Hole  features  are 
suppressed 


Round  features 
are  suppressed 


(b)  Suppression  of  features 


Figure  3.8:  Part  B:  boundary  condition  features  and  feature  suppression  after  rule 
implementation 
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Loading  is  applied 
to  the  entire  top  face 


conditions 


Round  features 
suppressed 


(b)  Suppression  of  features 


Figure  3.9:  Part  C:  boundary  condition  features  and  feature  suppression  after  rule 
implement  at  ion 


58 


(a)  Boundary  conditions 


(b)  Suppression  of  features 


Figure  3.10:  Part  D:  boundary  condition  features  and  feature  suppression  after  rule 
implement  at  ion 
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(a)  Boundary  conditions 


(b)  Suppression  of  features 


Figure  3.11:  Part  E:  boundary  condition  features  and  feature  suppression  after  rule 


implement  at  ion 
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Now,  an  analysis  was  conducted  independently  for  both  the  original  and  sim- 
plihed  models  of  these  parts  after  applying  the  boundary  conditions  described  above. 
First  the  learned  rules  were  applied  to  the  hve  test  parts.  After  applying  the  learned 
rules,  the  suppressibility .state  of  the  features  by  the  learned  rules  were  compared 
to  the  suppressibility .states  of  the  features  given  by  an  expert  analyst.  The  learned 
rules  proved  to  be  highly  accurate  for  the  hve  test  parts  compared  to  the  expert 
analysis  of  the  parts.  For  the  hve  parts  there  were  a  total  of  two  feature  classihcation 
errors  of  the  total  154  features  in  the  models,  both  arising  within  the  round  features. 
The  errors  were  caused  due  to  rare  suppression  occurrences  within  the  simplihcation 
method,  but  found  to  be  conservative  in  nature.  One  round  was  relatively  larger 
and  the  other  was  relatively  closer  to  a  boundary  condition  when  compared  to  the 
standard  values  for  round  suppression.  Table  3.7  displays  the  results  of  the  rules  to 
the  actual  suppression  analysis. 

Meshes  were  then  generated  for  these  parts  in  Pro/Mechanica  to  allow  for  a 
comparison  of  analysis  results  in  terms  of  mesh  quality,  mesh  generation  time,  and 
maximum  Von  Mises  stress  between  the  original  and  simplihed  models.  The  results 
shown  in  Table  3.8  demonstrate  a  signihcant  reduction  in  the  amount  of  elements 
by  more  than  a  factor  of  two,  with  the  largest  being  a  factor  of  over  eight.  There 
was  also  a  decrease  in  the  maximum  aspect  ratio,  signifying  a  higher  quality  mesh. 
Due  to  the  simplihcation,  the  time  to  mesh  was  also  greatly  reduced  by  at  least  27 
seconds;  this  could  be  signihcant  when  analyzing  numerous  parts. 

To  evaluate  the  ehectiveness  of  the  simplihcation  process,  the  maximum  Von 
Mises  stresses  were  compared  between  the  original  and  simplihed  models.  In  each 
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Feature  Type 

No.  of  Features 

No.  of  Errors 

Error  Type 

Hole 

93 

0 

None 

Round 

44 

2 

Conservative 

Chamfer 

17 

0 

None 

Table  3.7:  Analysis  of  learned  rules  on  the  five  test  parts  (A)-(E). 


case  there  is  a  difference  of  less  than  five  percent,  indicating  that  the  simplihcation 
process  did  not  significantly  alter  the  accuracy  of  the  results.  The  stress  analyses 
results  for  both  the  unsimplihed  and  simplihed  models  can  be  seen  in  Figs.  3.13- 
3.16.  In  the  results  of  part  C,  the  critical  stress  values  are  displayed  to  show  that 
this  location  is  constant  between  the  original  and  simplihed  models. 
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Tetra 

Max. 

Meshing  Time 

Max.  Von 

Elements 

Aspect  Ratio 

(sec) 

Mises  Stress  (MPa) 

Part-A  original 

19500 

11.39 

81 

33.97 

Part-A  simplified 

4511 

11.20 

7.2 

33.15 

Part-B  original 

6848 

20.09 

22.2 

372.39 

Part-B  simplified 

2462 

11.20 

7.8 

373.01 

Part-C  original 

20643 

20.34 

165 

30.06 

Part-C  simplified 

7865 

15.60 

46.9 

29.93 

Part-D  original 

7706 

11.17 

28.2 

42.27 

Part-D  simplified 

1445 

10.59 

4.8 

43.61 

Part-E  original 

9164 

11.45 

43.2 

162.30 

Part-E  simplified 

4049 

11.19 

13.8 

159.82 

Table  3.8:  Analysis  of  models  before  and  after  simplification 
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4 . 500e+03 
4 . 0006 ■•'03 
3.5006+03 
3 . 000e+03 
2.500e+03 
2 . 0006+03 
1 . 500e+03 
1 . 000e+03 
5 . 0006+02 


(a)  Unsimplified  model 


4 . 5006+03 
4 . 000e+03 
3.5006+03 
3 . 0006+03 
2.5006+03 
2. 000e+03 
1.5006+03 
1 . 0006+03 
5. 0006+02 

(PSI) 


(b)  Simplified  model 


Figure  3.12:  Stress  analysis  results  for  the  unsimplified  and  simplified  models  for 


part  A 
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4.86le+04 
4.320e+04 
3.7B0e+04 
3.2-40C  +  04 
2.700e+04 
2. 160e+04 
1.620e+04 
1 . 080e+0^ 
5 . 401e+03 


(PSI) 


(a)  Unsimplified  model 


(b)  Simplified  model 

Figure  3.13:  Stress  analyses  results  for  the  unsimplified  and  simplified  models  for 


part  B 
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4.360e+03 
3.924e+03 
3. 488e+03 
3.052e+03 
2.616e+03 
2. 100e+03 
1 . 744ff+03 
1 . 308e+03 
8.720i?  +  02 
A . 360e+02 
6.590e-02 
(PSI) 


(a)  Unsimplified  model 


4.342e+03 
3.924e+03 
3. 488e  +  03 
3.052e+03 
2.616e+03 
2. 100e+03 
1 . 744e+03 
1 .308e  +  03 
8.720e+02 
4. 360e+02 
1.792e-01 
(PSI) 


(b)  Simplified  model 


Figure  3.14:  Stress  analyses  results  for  the  unsimplified  and  simplified  models  for 
part  C 


66 


5.518e+03 
4.905e+03 
A . 291e+03 
3.678e+03 
3.065e+03 
2. 452e+03 
1 .839e+03 


(a)  Unsimplified  model 


5.694e+03 
5.062e+03 
4.42Se+03 
3.797e+03 
3. 165e+03 
2.533e+03 
1 . 900e+03 
1 .2680+03 
6. 360e+02 

(PSI) 


(b)  Simplified  model 

Figure  3.15:  Stress  analyses  results  for  the  unsimplified  and  simplified  models  for 
part  D 
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2. 119e+04 
1 .883e+04 
1 .648e+04 
1.413e+04 
1 . 177e+04 
9.417e+03 
7.063e+03 
4.709e+03 
2.354e+03 


a)  Unsimplified  model 


2.086e+04 
1.8540+04 
1 .622e+04 
1 .391e+04 
1 . 159e+04 
9.271e+03 
6.9530+03 
4.6360+03 
2.318e+03 


(b)  Simplified  model 


Figure  3.16:  Stress  analyses  results  for  the  unsimplified  and  simplified  models  for 


part  E 
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After  conducting  finite  element  analysis,  the  five  simplified  parts  showed  ac¬ 
ceptable  errors  within  5%  and  improved  mesh  quality  when  compared  with  the 
original  models.  The  learned  rules  proved  to  be  highly  accurate  for  the  test  parts. 
As  mentioned  earlier  only  conservative  errors  were  allowed  for  within  the  analysis 
process. 

Examination  of  how  the  parent-child  relationships  affect  the  results  was  done. 
For  that  the  parts  A  and  D  were  considered,  shown  in  Fig.  3.17(a)-(b).  Within 
these  models  the  parent  child  dependencies  were  accessed  and  there  were  features 
that  needed  to  be  retained  in  order  to  preserve  parent-child  relationships.  The 
chamfers  shown  in  red  in  Fig.  3.17(a)  has  its  children  shown  in  green.  If  the  parent 
is  suppressed  then  these  children  will  also  get  suppressed.  However,  in  this  case  the 
children  also  meet  the  suppression  criteria,  hence  both  the  parent  and  child  would 
be  suppressed.  The  holes  shown  in  the  Fig.  3.17(b)  cannot  be  suppressed  because 
the  children  are  labeled  BOUNDARY  ^CONDITION. 

(a)  Part-A  (b)  Part-D 

Figure  3.17:  Example  of  suppression  based  on  parent-child  relationships 
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3.5  Summary 


In  this  chapter,  a  framework  was  presented  that  can  be  utilized  in  any  of  the 
existing  CAD  softwares  for  suppressing  geometric  features  before  conducting  hnite 
element  analysis.  The  automated  tool  developed  in  this  chapter  assists  users  in 
identifying  the  suppressible  features  for  FEA  and  it  reduces  the  burden  on  users 
in  identifying  suppressible  features  manually.  Towards  generating  this  automated 
tool,  initially  a  rule  language  was  developed  that  can  capture  experts’  suppression 
knowledge.  This  knowledge  was  acquired  through  expert  demonstrations  and  the 
rules  were  expressed  in  terms  of  a  decision  tree.  With  the  help  of  statistical  induction 
learning  techniques  suppression  rules  were  learned  for  suppressing  holes,  rounds  and 
chamfers. 

The  Pro/Toolkit  API  available  in  ProE  was  utilized  to  generate  a  model  simpli¬ 
fication  tool;  for  hole,  round  and  chamfer  features.  The  complete  model  is  traversed 
through  to  determine  the  features  that  were  capable  of  suppression  based  upon  the 
rules  learned  through  expert  demonstrations.  The  suppression  state  was  determined 
by  the  feature  type,  dimensions,  and  relation  to  other  features,  i.e.  parent-child  re¬ 
lationship  and  inter-feature  distances.  To  demonstrate  the  simplihcation  method, 
five  different  models  were  analyzed.  Evaluation  of  the  learned  rules  found  two  out 
of  154  features  that  should  have  been  suppressed,  but  were  given  the  suppressibil- 
ity_state  of  FALSE;  this  error  was  of  conservative  nature.  In  comparison  to  the 
original  models,  the  simplihed  models  had  a  decrease  of  at  least  two  percent  in  the 
amount  of  elements  within  the  mesh.  Furthermore,  the  maximum  Von  Mises  stress 
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was  within  the  error  range  of  5%. 
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Chapter  4 


Assembly  Model  Simplification  for  Finite  Element  Analysis 
4.1  Overview 


The  main  focus  of  the  assembly  approach  is  to  design  a  framework  that  as¬ 
sists  the  user  in  the  simplification  of  complex  CAD  assembly  models,  consisting  of 
numerous  part  models;  similar  to  the  part  simplification  framework.  The  overall  ap¬ 
proach  can  be  broken  down  into  a  process  flow,  displayed  in  Fig.  4.1.  Initially,  the 
complete  assembly  information  is  extracted  and  written  in  an  XML  format.  This 
information  is  stored  as  an  element  tree  that  contains  all  information  about  the 
parts,  including,  for  example:  all  attributes,  such  as  the  material,  and  all  references 
to  the  existing  parts,  such  as  placement  and  orientation.  Upon  extraction  of  the 
model  information,  the  insignificant  parts  are  then  determined  and  the  model  is 
simplified,  whether  parts  are  suppressed  or  replaced. 


v4a  viIa 

Determtne  analysis  Verify  parent-child 

type  and  constraints  relationship  suppressions 


Figure  4.1:  Overview  of  the  assembly  approach 
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Discussions  with  AVMI  engineers  introduced  three  key  parameters  that  can  be 
used  to  assist  in  the  classihcation  of  parts  as  structurally  significant  or  insignificant. 
These  parameters  include  DESCRIP,  PART_NAME,  and  MATERIAL. 

1.  The  DESCRIP  parameter  often  provides  part  information.  This  information 
can  include  the  type  of  the  part;  e.g.  nut,  washer,  or  screw  (standard  hard¬ 
ware),  and  characteristics  of  the  part  type,  e.g.  thread  pitch,  bolt  diameter, 
and  head  shape  (bolt  or  screw). 

2.  The  generation  of  the  PART_NAME  is  usually  standardized  throughout  the 

industry,  grouping  similar  parts.  For  instance,  within  AVMI  a  PART_NAME 
including  the  AN*,  MS*,  or  NAS*  string  is  standard  hardware,  whereas  a 
PART_NAME  including  the  PE*  or  string  is  a  unique  part. 

3.  The  part  material  is  dependent  upon  the  requirements  of  the  part  within  the 
overall  assembly.  Different  materials  are  used  for  different  types  of  parts,  so  a 
specihc  material  can  be  seen  as  insignihcant  within  the  overall  model. 

A  search  within  the  XML  hie  for  key  words  and/or  numbering  schemes  of 
the  specihed  parameters  can  indicate  whether  or  not  the  part  is  standard  hardware 
or  a  unique  component.  Standard  hardware  is  generally  not  advantageous  within 
the  FEA  model,  as  they  generate  unnecessary  mesh  complexities  within  the  analysis 
and  can  be  suppressed;  whereas  a  unique  component  is  a  designed  part  that  requires 
FEA,  but  may  be  capable  of  simplihcation. 

The  framework  has  been  implemented  within  the  ProE  CAD  system.  Section 
4.2  describes  the  algorithm  for  suppressing  parts  based  upon  the  PART  .NAME  and 
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MATERIAL  and  usefulness  of  this  technique.  Section  4.3  describes  simplihcation 
methods  pertaining  to  the  individual  parts  of  the  assembly.  Section  4.4  presents 
results  of  the  simplihcation  process,  by  using  the  techniques  described  in  the  previous 
sections. 

4.2  Part  Suppression 

In  order  to  assist  in  the  simplihcation  process  of  the  assembly  models,  the 
Pro/Toolkit  API  in  ProE  was  utilized  to  generate  a  model  suppression  tool,  just 
as  in  the  part  simplihcation  framework.  For  recap,  this  tool  is  an  add-on  module 
in  ProE,  accessible  through  the  menu  systems  in  the  ProE  environment.  In  the 
API,  the  design  model  is  stored  by  ProMdl,  a  data  handle  to  retrieve  the  complete 
information  of  the  design  model.  This  information  is  stored  within  ProE  as  an 
element  tree  that  contains  all  information  about  the  parts,  including,  for  example: 
all  attributes,  such  as  the  material,  and  all  references  to  the  existing  parts,  such  as 
placement  and  orientation.  Suppression  is  based  upon  the  part  MATERIAL  and 
PART_NAME  parameters. 

4.2.1  Part  Utilization 

Dehning  the  intended  purpose  of  the  part  within  the  overall  assembly  model; 
e.g.  structural  stability,  thermal  enhancement,  ergonomic  aspect,  or  part  attach¬ 
ment  method;  can  facilitate  the  classihcation  of  the  part’s  criticality  for  analysis. 
The  PART_NAME  and  DESCRIP  parameters  can  aid  in  the  identihcation  of  the 
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part’s  purpose  within  the  assembly.  Generally,  the  PART_NAME  has  a  common 
string,  specifying  the  group  type  to  which  the  part  belongs;  for  example  the  MS* 
string  refers  to  a  type  of  fastener  and  then  the  DESCRIP  parameter  contains  key 
words  to  assist  in  specifying  the  intended  purpose,  such  as  bolt  or  screw. 

Although  this  can  facilitate  the  classihcation,  the  intended  purpose  of  the 
part  is  not  always  clearly  dehned.  Different  aspects  of  the  part  can  be  intended  for 
various  purposes  or  parts  unintended  for  the  specihed  purpose  can  have  an  effect 
on  different  analysis  results;  e.g.  a  heat  shield  that  is  placed  within  an  assembly 
for  the  thermal  protection  of  other  parts,  may  not  only  affect  a  thermal  analysis 
(for  which  the  part  is  intended),  but  could  also  have  an  effect  on  the  modal  analysis 
of  the  overall  model.  For  these  reasons,  focus  was  geared  toward  a  specific  part 
type  that  mainly  falls  into  the  intended  purpose  of  attachment  methods:  hardware 
(i.e.  fasteners),  such  as  nuts  and  bolts.  Depending  on  the  number  of  parts  in  the 
assembly,  a  model  could  potentially  contain  a  couple  hundred  fasteners.  Even  with 
the  large  amount  of  fasteners,  there  is  often  only  a  select  few  of  different  fasteners 
within  the  model.  Using  the  same  fasteners  throughout  reduces  the  overall  part 
count  of  the  model.  There  are  several  benefits  to  reducing  the  overall  part  count 
of  an  assembly;  (l)decreases  some  of  the  cost  because  less  material  types  need  to 
be  purchased,  (2)  similar  tools  can  be  used  to  assemble  the  model  and  (3)  reduces 
the  error  probability  during  the  manufacturing  process  of  the  model.  Typically, 
these  fasteners  have  minimal  to  no  effect  on  analyses  and,  therefore,  fasteners  are 
frequently  suppressed  before  an  analysis;  greatly  reducing  the  computational  cost 
of  the  analysis. 
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Currently,  an  engineer  will  have  to  manually  transverse  through  the  model 
and/or  the  model  tree  and  independently  suppress  each  fastener.  With  a  large 
model  this  could  become  a  tedious  time  consuming  task.  To  help  expedite  the 
process  of  fastener  suppression,  or  other  parts  that  may  appear  several  times  within 
the  model,  the  framework  can  suppress  the  parts  by  the  PART  .NAME;  all  parts  in 
the  model  with  that  part  name  will  be  quickly  suppressed.  For  instance,  a  model 
containing  25  bolts  of  two  different  types  (PART_NAMEs),  would  only  require  two 
independent  suppressions,  instead  of  the  manual  independent  suppressions  of  25. 

4.2.2  Material 

In  addition  to  the  part’s  intended  purpose,  material  selection  plays  a  role  in 
determining  the  extent  of  the  part  criticality,  from  the  analysis  results  point  of 
view.  The  analysis  response  of  a  part  is  highly  dependent  upon  its  overall  shape 
and  dimensions,  but  the  magnitude  of  this  response  is  a  result  of  the  material;  i.e. 
the  material  properties.  For  instance.  Fig.  4.2  displays  an  example  assembly  model 
section,  which  is  the  input  for  two  static  analyses.  For  comparison  the  parts  in  red 
are  constucted  of  a  similar  material  and  only  the  material  of  these  two  parts  is  altered 
between  the  two  tests.  The  material  properties  of  the  two  parts  for  the  analyses  are 
that  of  steel  and  nylon.  From  the  results  shown  in  Figs.  4.3  and  4.4,  steel  and  nylon 
material  attributes  respectively,  it  can  be  deduced  that  the  critical  stress  within 
the  model  is  much  higher  in  the  assembly  model  where  the  parts  are  made  of  steel, 
than  that  of  the  assembly  model  where  the  parts  are  made  of  nylon.  This  stress 
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Figure  4.2:  Example  assembly  model 
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Figure  4.3:  Stresses  within  the  model  when  the  two  parts  are  made  of  the  material 


steel 


difference  is  the  result  of  the  difference  in  the  densities  of  the  two  materials.  The 
analysis  demonstrates  that  the  part  material  is  critical  in  determining  the  criticality 
of  the  part  for  analysis  purposes. 

The  behavioral  response  of  the  material,  due  to  its  properties  as  discussed 
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above,  has  generated  basic  standards  within  the  material  selection  phase;  for  in¬ 
stance,  a  load  bearing  part  is  generally  constructed  of  a  material  with  a  higher 
strength,  such  as  a  metal  compared  to  a  plastic.  Knowledge  of  these  standards  and 
the  analysis  type  to  be  performed  on  the  model  can  signify  material  types  as  in- 
signihcant  or  signihcant  to  the  analysis.  For  example,  in  a  structural  analysis,  often, 
plastic  parts  can  be  viewed  as  insignihcant,  since  load  bearing  parts  are  usually  not 
constructed  with  this  material. 


4.2.3  Algorithm 

The  parts  considered  to  be  insignihcant  based  upon  their  material  or  purpose 
are  then  suppressed.  Either  the  MATERIAL  or  PART_NAME  value  can  be  entered 
to  suppress  all  the  parts  containing  that  parameter  value.  Before  suppression  it 
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is  necessary  to  understand  the  dependency  relationships  among  the  parts,  just  as 
discussed  in  the  previous  chapter.  The  previous  chapter  discusses  the  feature-based 
approach  in  model  generation,  where  features  are  dependent  among  one  another 
within  a  model;  generating  the  feature  history  tree  [33,  37].  Just  as  the  features  are 
dependent  upon  one  another  in  a  part  model,  parts  are  dependent  upon  one  another 
within  an  assembly  model.  When  constructing  the  assembly  model,  parts  are  intro¬ 
duced  into  the  model  one  at  a  time  referencing  the  previous  parts  within  the  model 
for  placement.  The  same  model  can  be  constructed  in  several  ways  ,  dependent 
upon  the  order  in  which  the  parts  were  added  to  the  model,  generating  different 
history  trees  and,  ultimately,  different  dependencies  among  the  parts.  Suppression 
of  a  parent  will  ultimately  suppress  its  children  parts,  hence,  it  is  imperative  to 
understand  parts’  parent-child  relationships  for  suppressions;  since  the  child  could 
potentially  be  a  critical  part  for  the  analysis. 

For  example,  let’s  look  at  the  construction  of  a  simple  table  with  four  legs, 
attached  by  nuts  and  bolts,  illustrated  in  Fig.  4.5.  Conventionally  the  table  top 
would  be  first  inserted  into  the  model.  Then  the  legs  would  be  brought  in  one  at  a 
time,  referencing  the  axis  of  the  holes  within  the  table  top  for  the  desired  location  of 
the  legs.  Upon  including  the  legs,  the  bolts  would  be  inserted  into  the  table  top  and 
legs  by  referencing  the  axes  once  again  to  line  up  the  bolts  and  then  referencing  the 
table  top  surface  and  connecting  surface  of  the  bolt  head.  After  this  the  nuts  are 
added  by  referencing  the  bolt  axis  and  the  connecting  surfaces  of  the  nut  and  table 
leg.  In  this  instance,  suppression  of  the  bolts  would  also  suppress  the  nuts,  since  they 
are  referenced  from  the  bolts  axis.  This  setup  would  not  cause  any  critical  features 
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Figure  4.5:  Simple  four-legged  table  for  illustration  of  the  parent-child  relationship 
within  a  model 

to  be  suppressed  within  the  model,  however,  a  designer  could  assemble  the  table 
in  a  different  manner  within  the  software.  Let’s  say  the  designer  unconventionally 
decides  to  bring  the  bolts  into  the  model  and  reference  the  table  top  and  legs  from 
the  bolts.  Then  suppression  of  the  bolts  would  ultimately  suppress  the  entire  model. 
Extracting  these  parent-child  dependencies  is  crucial  in  order  to  avoid  suppressing 
any  non-critical  features. 

An  algorithm  has  been  developed  to  account  for  the  parent-child  relationships 
during  the  suppression  process,  similar  to  that  of  the  algorithm  presented  in  the 
previous  chapter  for  suppressing  features  based  upon  the  parent-child  relationships. 
Alteration  of  the  algorithm  was  required,  since  there  is  not  a  general  rule  set  that 
assigns  a  suppressibility  status  based  upon  all  the  criteria;  i.e.  suppression  is  based 
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upon  either  the  MATERIAL  or  PART_NAME  value. 

The  algorithm  for  suppressing  features  is  as  follows: 

1.  Construct  G  (parent-child  relationship  graph)  by  inserting  a  node  V  corre¬ 
sponding  to  every  part  in  the  tree,  which  is  not  suppressed.  If  part  p’  uses 
a  part  p  as  a  reference,  then  insert  an  edge  (v,  p’)  in  G.  Here  the  vertex  v 
corresponds  to  part  p  and  vertex  v’  corresponds  to  part  p’.  Fig.  4.6  shows  an 
example  of  the  parent-child  relationship  graph. 

2.  Find  nodes  X  in  V  that  do  not  have  children  and  have  the  attribnte  valne  (MA¬ 
TERIAL  or  PART  .NAME)  given.  If  X  is  empty  proceed  to  step  3,  otherwise 
do  the  following: 

(a)  Remove  every  node  in  X  from  G  by  deleting  the  associated  vertex  and 
edges. 

(b)  Go  to  step  2. 

3.  Find  nodes  Y  in  V  that  have  children  and  the  attribnte  value  and  do  the 
following: 

(a)  For  node  Yi,  highlight  the  child/children  in  green  and  the  parent  node 
in  red,  ask  if  the  user  would  like  to  suppress  the  parent  and  associated 
child/children.  If  Yi  is  empty,  then  stop. 

i.  Yes,  then  remove  node  Yi  and  child/children  nodes  of  G  by  deleting 
the  associated  vertex  and  edges.  Proceed  back  to  step  3a. 
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Figure  4.6:  Example  of  the  generated  parent-child  relationship  graph 


ii.  No,  then  add  1  to  i  and  proceed  to  step  3a. 

This  model  is  only  valid  for  assembly  models  generated  using  only  parts  and  no 
minor  assembly  models.  When  minor  assembly  models  are  used  within  the  overall 
assembly,  the  parent  references  of  the  minor  assembly  are  not  necessarily  used  as 
the  reference  for  inclusion  into  the  overall  model.  If  the  parent  is  not  used,  then 
it  is  not  seen  as  a  child  to  the  parts  from  which  the  minor  assembly  is  referenced. 
Then  suppression  of  the  part  model  of  the  overall  assembly  used  for  referencing 
generates  an  error  because  the  parents  of  the  referenced  part  are  not  suppressed, 
but  the  references  for  the  part  into  the  overall  model  are  suppressed. 

4.3  Simplification 

Suppression  of  parts  is  not  the  only  means  of  simplifying  an  assembly  model, 
the  parts  themselves  can  also  be  simplified  within  the  model.  Simplification  of  the 
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parts  can  be  accomplished  by  suppressing  the  non-critical  features,  or  the  parts  can 
be  replaced,  either  with  a  simpler  geometry  or  point  mass  that  maintains  the  effects 
of  the  original  part. 

The  hrst  method  was  thoroughly  discussed  in  Chapter  2.  To  successfully  use 
this  approach  within  ProE,  the  simplihcation  of  the  part  must  be  done  within  the 
part  hie  because  modihcations  cannot  be  done  to  an  individual  part  within  the 
assembly  hie.  Even  though  the  modihcations,  such  as  suppression  of  features,  are 
done  within  the  part  hie,  they  will  be  automatically  rehected  within  the  assembly 
model. 

The  second  method  requires  additional  knowledge  compared  to  that  of  the 
other  method,  since  the  features  are  not  only  suppressed,  but  the  part  geometry 
is  altered;  maintaining  a  similar  response  to  that  of  the  original  part.  First,  the 
analysis  type  must  be  identihed,  since  a  part  can  have  a  diherent  response  within 
the  various  types  of  analyses.  Upon  knowledge  of  the  analysis  type,  the  part  must 
be  quickly  reviewed  to  determine  the  necessity  and  response,  so  a  suitable  approx¬ 
imation  method  can  be  developed,  whether  the  method  is  point  mass  or  simpler 
geometry  with  density  match.  Then  the  center  of  gravity  (CG)  is  calculated  of  the 
original  part.  For  the  point  mass  simplihcation  the  CG  is  the  location  of  the  point 
mass.  For  the  simpler  geometry  simplihcation  the  simpler  geometry  is  designed  to 
maintain  the  location  of  the  CG.  Figure  X  illustrates  a  part  that  was  simplihed 
to  a  simpler  geometry  displayed  in  Figure  Y.  To  automate  this  method  an  outside 
program  requiring  numerous  inputs  will  be  needed  to  analyze  and  create  a  simpler 
geometry  according  to  the  model. 
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Figure  4.7:  Unsimplified  boom  model  for  a  modal  analysis 

4.4  Results  and  Discussion 

To  illustrate  the  framework  described,  assembly  models  have  been  simplified 
for  both  a  structural  and  a  modal  analysis.  Simplihcation  of  the  models  was  accom¬ 
plished  by  part  suppression,  based  upon  the  MATERIAL  or  PART_NAME  attribute 
value,  and  part  simplihcation  discussed  in  the  previous  chapter.  Upon  model  sim- 
plihcation  a  hnite  element  analysis  is  conducted  for  both  the  simplihed  and  original 
models,  with  the  help  of  Pro/Mechanica.  After  conducting  the  FEA,  the  analysis 
results  were  compared  between  the  simplihed  and  original  models  in  terms  of  mesh 
quality  (number  of  tetra  elements  and  maximum  aspect  ratio),  mesh  generation  time, 
and  either  mode  values  or  maximum  Von  Mises  stress,  depending  on  the  analysis 
type. 

4.4.1  Modal  Analysis 

A  modal  analysis  was  done  on  a  boom  assembly  model  shown  in  Fig.  4.7.  The 
simplihcation  of  the  assembly  model,  displayed  in  Fig.  4.8,  was  done  in  the  following 
manner: 

•  The  15  bolts  and  their  respective  nuts  were  suppressed  by  their  PART_NAME. 
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Figure  4.8:  Simplified  boom  model  for  a  modal  analysis 


To  suppress  the  15  nuts  it  only  took  two  different  PART_NAMEs,  which  ac¬ 
counted  for  the  two  different  diameter  nuts;  although  it  required  four  PART_NAMEs 
to  suppress  the  bolts,  accounting  for  the  different  lengths  and  diameters.  To 
complete  the  suppression  of  the  30  parts  it  only  took  6  passes  compared  to  the 
manual  suppression  of  30  passes.  For  clarihcation  the  nuts  were  suppressed 
prior  to  the  bolts,  removing  the  children  parts  before  the  parent  parts.  If  the 
bolts  were  hrst  suppressed,  then  the  nuts  would  have  been  seen  as  children  to 
the  bolts  and  the  suppression  message  would  have  appeared.  Either  order  of 
suppression  achieves  the  same  desired  outcome. 

•  The  two  mount  attachments,  the  center  slider,  the  complex  attachment,  and 
the  boom  were  simplihed  using  the  method  discussed  in  the  previous  chapter. 

The  hrst  step  in  the  FEA  process  is  mesh  generation.  The  comparison  of  the 
mesh  quality  and  mesh  generation  time  between  the  two  models  can  be  seen  in  Table 
4.1.  As  illustrated  below  the  simplihed  model  has  a  44%  reduction  of  the  number 
of  elements;  the  max  aspect  ratio  is  reduced  by  7.75,  which  means  the  mesh  quality 
is  better  in  the  simplihed  model;  and  the  time  to  mesh  was  reduced  by  almost  20% 
or  13.8  sec.  In  terms  of  the  mesh,  the  simplihed  model  has  a  better  quality  mesh 
and  was  generated  faster  than  the  original  model. 
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Tetra  Elements 

Max.  Aspect  Ratio 

Meshing  Time  (sec) 

Original 

15069 

20.90 

72 

Simplified 

8422 

13.15 

58.2 

Table  4.1;  Mesh  improvement 


Mode 

Original  Model  (Hz) 

Simplified  Model  (Hz) 

%  Error 

1 

25.23 

25.28 

0.1982 

2 

73.50 

73.47 

0.0408 

3 

75.32 

74.75 

0.7568 

4 

77.21 

77.16 

0.0648 

Table  4.2:  Modes  calculated  within  the  modal  analysis 


The  mesh  quality  is  not  the  only  factor  required  for  comparison  of  the  models, 
but  the  analysis  results  of  the  model  must  be  similar  or  the  simplihed  model  cannot 
replace  the  original  model  within  the  analysis  process.  A  modal  analysis  has  been 
executed  on  the  boom  models  and  the  hrst  four  modes  of  both  the  original  and 
simplihed  models  were  compared;  shown  in  Table  4.2.  The  last  column  of  the  table 
shows  that  there  was  an  error  of  less  than  1%  between  the  two  models  for  all  four 
of  the  modes. 


The  modes  of  the  two  models  can  be  seen  in  Figs.  4.9-4.12. 
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Figure  4.9:  (a)  Mode  1  of  the  unsimplified  boom  assembly  (b)  Mode  1  of  the  sim¬ 


plified  boom  assembly 
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Figure  4.10:  (a)  Mode  2  of  the  unsimplified  boom  assembly  (b)  Mode  2  of  the 


simplified  boom  assembly 
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Figure  4.11:  (a)  Mode  3  of  the  unsimplified  boom  assembly  (b)  Mode  3  of  the 


simplified  boom  assembly 
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Figure  4.12:  (a)  Mode  4  of  the  unsimplified  boom  assembly  (b)  Mode  4  of  the 
simplihed  boom  assembly 
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4.4.2  Structural  Analysis 


In  a  similar  fashion  to  the  modal  analysis,  a  structural  analysis  was  completed 
for  two  assembly  models.  The  two  models  used  for  the  structural  analysis  were  a 
weight  bench  assembly  and  the  same  boom  assembly  from  which  the  modal  analysis 
was  done.  Comparison  between  the  simplihcation  process  for  the  modal  analysis 
and  structural  analysis  of  the  beam  will  highlight  different  simplihcation  processes. 

First,  the  original  and  simplihed  models  of  the  weight  bench  are  shown  in 
Figs.  4.13  and  4.14,  respectively.  Simplihcation  of  the  assembly  was  done  in  the 
following  manner: 


Figure  4.13:  Unsimplihed  weight  bench  assembly  model  for  a  structural  analysis 


•  The  12  bolts  and  their  respective  nuts  were  suppressed  by  their  PART_NAME. 
To  suppress  the  12  nuts  it  only  took  two  diherent  PART_NAMEs,  which 
accounted  for  the  two  diherent  diameter  nuts;  although  it  required  three 
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Figure  4.14:  Simplified  weight  bench  assembly  model  for  a  structural  analysis 

PART_NAMEs  to  suppress  the  bolts,  accounting  for  the  different  lengths  and 
diameters.  To  complete  the  suppression  of  the  24  parts  it  only  took  five  passes 
compared  to  the  manual  suppression  of  24  passes.  For  clarification  the  nuts 
were  suppressed  prior  to  the  bolts,  removing  the  children  parts  before  the 
parent  parts.  If  the  bolts  were  first  suppressed,  then  the  nuts  would  have 
been  seen  as  children  to  the  bolts  and  the  suppression  message  would  have 
appeared.  Either  order  of  suppression  achieves  the  same  desired  outcome. 

•  The  support  frames  were  simplified  using  the  method  discussed  in  the  previous 
chapter. 

•  The  attachments  covering  the  ends  of  the  pipes  are  all  constructed  of  a  nylon 
material  and  have  no  bearing  on  the  stresses  generated  within  the  model. 
These  9  caps  were  suppressed  in  one  pass,  using  their  MATERIAL  value. 
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Just  as  in  the  modal  analysis,  the  first  step  in  the  structural  analysis  is  mesh 
generation.  The  comparison  of  the  mesh  quality  and  mesh  generation  time  between 
the  two  models  can  be  seen  in  Table  4.3.  As  illustrated  below  the  simplihed  model 
has  a  35%  reduction  of  the  number  of  elements;  the  max  aspect  ratio  is  reduced  by 
2.27,  which  means  the  mesh  quality  is  better  in  the  simplihed  model;  amd  the  time  to 
mesh  was  rednced  by  abont  38%  or  150.68  sec.  In  terms  of  the  mesh,  the  simplihed 
model  has  a  better  qnality  mesh  and  was  generated  faster  than  the  original  model. 
Upon  mesh  generation  the  Von  Mises  stresses  are  calculated  within  the  model  and 
the  maximnm  stress  is  then  compared  between  the  two  models,  shown  in  Table  4.3, 
which  has  an  error  of  4.8%  between  the  original  and  simplihed  weight  bench  models. 
Figs.  4.15  and  4.16  show  the  stresses  within  the  two  models 


Tetra 

M2Lximum 

Meshing  Time 

Max.  Von 

Elements 

Aspect  Ratio 

(sec) 

Mises  Stress 

Original 

46423 

19.78 

391.8 

39.93 

Simplified 

30133 

17.51 

241.2 

41.85 

Table  4.3:  Mesh  improvement 
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Figure  4.15:  Von  Mises  Stress  distribution  on  the  unsimplified  weight  bench  assem¬ 


bly 


Figure  4.16:  Von  Mises  Stress  distribution  on  the  simplified  weight  bench  assembly 


94 


The  other  assembly  model  for  which  a  sructural  analysis  was  done  was  the 
boom  assembly.  Figs.  4.17  and  4.18  illustrate  the  original  and  simplihed  models 
of  the  boom  assembly.  The  original  model,  Fig.  4.17,  is  the  same  as  the  original 
modal  analysis  model.  Fig.  4.7,  but  the  simplihed  models  are  slightly  altered.  The 
structural  analysis  simplihed  model  has  replaced  the  complex  structure  on  the  end 
of  the  boom  with  a  simpler  part;  where  in  the  modal  analysis  this  complex  end  was 
not  replaced,  but  just  simplihed,  according  to  the  previous  chapter.  Simplihcation 
of  the  assembly  was  done  in  the  following  manner: 

Figure  4.17:  Unsimplihed  boom  assembly  model  for  a  structural  analysis 


Figure  4.18:  Simplihed  boom  assembly  model  for  a  structural  analysis 


•  The  15  bolts  and  their  respective  nuts  were  suppressed  by  their  PART_NAME. 

To  suppress  the  15  nuts  it  only  took  two  diherent  PART_NAMEs,  which  ac¬ 
counted  for  the  two  diherent  diameter  nuts;  although  it  required  four  PART_NAMEs 
to  suppress  the  bolts,  accounting  for  the  diherent  lengths  and  diameters.  To 
complete  the  suppression  of  the  30  parts  it  only  took  6  passes  compared  to  the 
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manual  suppression  of  30  passes.  For  clarification  the  nuts  were  suppressed 
prior  to  the  bolts,  removing  the  children  parts  before  the  parent  parts.  If  the 
bolts  were  hrst  suppressed,  then  the  nuts  would  have  been  seen  as  children  to 
the  bolts  and  the  suppression  message  would  have  appeared.  Either  order  of 
suppression  achieves  the  same  desired  outcome. 

•  The  two  mount  attachments,  the  center  slider,  and  the  boom  were  simplihed 
using  the  method  discussed  in  the  previous  chapter. 

•  The  complex  attachment  at  the  end  was  replaced  by  a  simpler  part  geometry. 
The  complexity  of  the  part  was  reduced  to  a  cylindrical  geometry,  maintaining 
the  CG  of  the  part.  This  was  not  capable  within  the  modal  analysis  because 
the  as  the  mode  number  increased  there  was  a  higher  emphasis  on  the  mode 
value  from  the  complex  features  of  the  attachment. 

The  comparison  of  the  mesh  quality  and  mesh  generation  time  between  the 
two  models  can  be  seen  in  Table  4.4.  As  illustrated  below  the  simplihed  model  has 
a  35%  reduction  of  the  number  of  elements;  the  max  aspect  ratio  is  reduced  by  9.23, 
which  means  the  mesh  quality  is  better  in  the  simplihed  model;  and  the  time  to 
mesh  was  reduced  by  about  69.2%  or  49.8  sec.  In  terms  of  the  mesh,  the  simplihed 
model  has  a  better  quality  mesh  and  was  generated  faster  than  the  original  model. 
Upon  mesh  generation  the  Von  Mises  stresses  are  calculated  within  the  model  and 
the  maximum  stress  is  then  compared  between  the  two  models,  shown  in  Table  4.4, 
which  has  an  error  of  4.8%  between  the  original  and  simplihed  boom  models.  Figs. 
4.19  and  4.20  show  the  stresses  within  the  two  models. 
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Tetra 

M^Lximum 

Meshing  Time 

Max.  Von 

Elements 

Aspect  Ratio 

(sec) 

Mises  Stress 

Original 

15069 

20.90 

72 

101.49 

Simplified 

5141 

11.67 

22.2 

96.59 

Table  4.4:  Mesh  improvement 


Figure  4.19:  Von  Mises  Stress  distribution  on  the  unsimplihed  boom  assembly 


Figure  4.20:  Von  Mises  Stress  distribution  on  the  unsimplihed  boom  assembly 
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4.5  Summary 


This  chapter  has  presented  a  tool  to  be  used  for  the  simplification  of  an  assem¬ 
bly  model  for  FEA,  utilizing  the  Pro/Toolkit  available  within  ProE.  The  user  must 
manually  determine  the  criticality  of  the  parts  within  the  model,  whether  they  are 
suppressible  or  non-suppressible  and  whether  the  non-suppressible  features  are  ca¬ 
pable  of  simplihcation,  based  upon  the  analysis  type,  material,  and  the  part  feature. 
Simplihcation  of  the  part  can  either  be  accomplished  by  the  method  proposed  in  the 
previous  chapter  or  by  replacing  the  complex  part  with  a  simpler  geometrical  part. 
After  the  criticality  of  the  part  is  determined,  the  suppressible  parts  can  be  sup¬ 
pressed  in  groups,  rather  than  one-by-one,  based  upon  either  their  PART_NAME  or 
MATERIAL  value.  This  also  takes  into  account  the  parent-child  relationships  devel¬ 
oped  within  the  model.  To  demonstrate  the  process,  two  different  assembly  models 
were  simplihed  and  analyzed.  In  comparison  to  manually  simplifying  the  parts,  the 
number  of  passes  to  suppress  the  non-critical  parts  was  better  than  halved;  while 
still  increasing  the  mesh  quality,  reducing  the  mesh  generation  time,  and  acquiring 
analysis  results  within  a  given  error  of  5%. 


Chapter  5 


Conclusion 

5.1  Intellectual  Contributions 

The  contributions  of  this  work  can  be  separated  into  two  primary  categories: 

1.  A  framework  for  simplifying  CAD  part  models  was  developed  with  the  utili¬ 
ties  provided  within  ProE.  The  approach  contains  several  key  aspects.  First, 
numerous  features  can  be  suppressed  at  one  time,  compared  to  the  industry’s 
current  simplihcation  process  of  suppressing  one  feature  at  a  time.  Second,  a 
rule  basis  has  been  learned  through  expert  demonstration  with  only  conserva¬ 
tive  error,  which  can  evolve  as  more  parts  are  simplihed.  Third,  the  parent- 
child  relationships  were  extracted  from  the  model  and  used  in  the  analysis  to 
determine  the  suppressibility  state  of  the  model,  thus,  ensuring  that  a  child 
feature  was  not  unintentionally  suppressed  within  the  process.  Finally,  in  or¬ 
der  to  reduce  the  potential  effect  of  the  simplihcation  on  the  FEA  results,  the 
inter-feature  distance  to  the  boundary  conditions  were  calculated  to  maintain 
any  features  within  a  close  proximity  of  the  boundary  and  force  constraints. 

2.  A  framework  for  simplifying  CAD  assembly  models  was  developed  with  the 
utilities  provided  within  ProE.  The  approach  contains  several  key  aspects. 
First,  numerous  parts  can  be  suppressed  at  one  time,  compared  to  the  indus- 
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try’s  current  simplification  process  of  suppressing  one  part  at  a  time.  Second, 
the  parent-child  relationships  were  extracted  from  the  model  and  used  when 
suppressing  the  parts,  ensuring  that  a  child  part  was  not  unintentionally  sup¬ 
pressed  within  the  process.  Finally,  integration  of  the  part  simplihcation  was 
introduced  within  the  assembly  model  simplihcation  to  simplify  parts  within 
the  overall  assembly  model. 

The  work  presented  in  this  thesis  provides  a  simplihcation  tool  for  both  part 
and  assembly  models,  which  is  not  the  case  for  all  model  simplihcation  tools. 

5.2  Anticipated  Benefits 

The  simplihcation  frameworks  can  be  extremely  benehcial  within  the  industry 
analysis.  The  overall  analysis  time  for  a  model  will  likely  be  reduced,  since  the  model 
preparation  time  can  be  reduced.  In  the  part  simplihcation  approach  discussed 
in  Chapter  2,  the  process  greatly  reduces  the  users  involvement  in  the  process. 
Automatic  identihcation  and  labeling  of  the  non-critical  features  and  multi-feature 
suppression  of  the  identihed  non-critical  features,  upon  the  user’s  verihcation,  will 
show  a  signihcant  time  reduction  from  the  manual  technique  of  the  user  analyzing 
each  feature  and  suppressing  them  one  at  a  time.  In  addition  to  reducing  the  model 
preparation  time,  errors  involved  in  the  simplifation  part  models  may  be  reduced; 
especially  in  the  accidental  suppression  of  critical  child  features,  since  the  the  parent- 
child  relationships  are  incorporated  into  the  suppression  algorithm.  Reducing  the 
errors  within  the  analyses  has  two  benehts.  First  it  reduces  the  chance  of  passing 
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along  erroneous  results  and  the  need  to  rerun  an  analysis  because  it  was  discovered 
to  have  suppressed  a  critical  part. 

In  the  assembly  model  simplihcation  tool,  numerous  parts  are  capable  of  sup¬ 
pression  compared  to  the  industry  version  of  suppressing  part  by  part;  reducing  the 
model  preparation  time.  Also,  the  integration  of  the  part  simplihcation  tool  into 
the  assembly  model  simplihcation  will  reduce  the  time  to  simplify  the  parts  within 
the  overall  model.  Initially,  these  tools  will  not  create  a  large  impact  in  reducing 
the  analysis  time,  but  as  time  goes  on  there  will  be  a  gradual  decrease  in  the  model 
preparation  time;  since  it  will  take  time  to  gain  the  engineers  trust. 

5.3  Future  Work 

While  this  thesis  has  described  successful  frameworks  for  the  simplihcation  of 
both  part  and  assembly  CAD  models  for  the  use  of  downstream  applications,  such 
as  FEA,  there  are  many  areas  in  which  future  work  can  oher  further  insight.  The 
main  areas  of  interest  for  future  improvement  are  directly  related  to  the  work  that 
has  been  complete  in  this  thesis. 

5.3.1  Improvement  to  the  Part  Simplification  Framework 

The  part  simplihcation  framework  presented  in  Chapter  2,  but  has  several  lim¬ 
itations.  First,  the  framework  is  limited  to  suppressing  only  three  types  of  features 
within  the  model;  holes,  rounds,  and  chamfers.  Within  the  held,  simplihcation  of  a 
part  is  not  only  limited  to  these  three  features.  Expansion  of  the  framework  to  in- 
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elude  other  features,  for  instance,  bosses,  slots,  and  steps,  would  greatly  improve  the 
robustness  of  the  overall  simplification  tool.  Another  aspect  of  the  framework  that 
could  be  improved  is  the  calculation  of  the  feature  distance  to  the  boundary  con¬ 
ditions.  In  the  current  framework  the  inter-feature  distance  calculation  only  takes 
into  account  the  geometric  distance.  Within  an  analysis  the  topological  distance 
and  stress  path  can  also  have  an  effect  within  the  criticality  of  the  feature,  so  these 
should  be  incorporated  into  or  included  with  the  geometric  distance  calcultaion. 

A  third  aspect  for  expansion  is  to  incorporate  the  overall  volnmetric  change 
between  the  original  and  simplihed  models.  The  overall  volnmetric  change  is  sim¬ 
ilar  to  the  parent-child  relationship,  in  the  fact  that  it  may  reqnire  features  with 
parameters  qualifying  it  to  be  a  snppressible  featnre  to  be  maintained  within  the 
model.  For  instance,  one  hole  having  parameters,  which  constitute  the  feature  to  be 
cionsidered  non-critical,  will  not  have  an  impact  on  the  analysis  results;  although, 
taking  this  same  hole  and  placing  it  nnmerous  times  throughont  the  model,  gen¬ 
erating  a  large  volumetric  change  between  the  simplihed  and  non-simplihed  model 
may  affect  the  analysis  results  (even  though  the  parameters  of  the  holes  constituted 
it  to  be  classihed  as  a  non-critical  featnre). 

5.3.2  Advancements  in  Assembly  Simplification  Approach 

The  assembly  model  simplihcation  framework  was  more  user  intensive  than 
that  of  the  part  simplihcation  framework,  allowing  for  more  areas  of  work  to  be 
focused  on  in  the  future.  First,  the  parent-child  relationship  would  need  to  be 
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extended  to  allow  for  the  inclusion  of  minor  assembly  models  within  the  overall 
assembly  model.  If  the  parent  part  of  an  assembly  model  to  be  included  is  not  used 
as  the  references  for  placement  into  the  model,  then  they  are  technically  a  child  in 
the  major  assembly,  for  which  they  are  a  parent  to  in  the  minor  assembly  model. 
Although  this  is  the  case,  the  possibility  exists  for  minor  assembly  parts  to  become 
unreferenced  in  the  overall  model,  if  its  child  is  used  as  the  placement  reference  for 
the  minor  assembly  and  the  suppressed  for  simplihcation  in  the  overall  model. 

Another  possible  area  for  extension  of  the  framework  is  to  automatically  de¬ 
termine  the  criticality;  non-critical,  semi-critical,  or  critical;  since  in  the  current 
framework  the  user  must  manually  determine  the  criticality  of  each  part.  Deter¬ 
mining  the  criticality  of  parts  requires  a  knowledge  of  not  only  the  part  attributes, 
but  also  of  the  assembly  attributes  and  analysis  type.  Chapter  4  illustrated  that 
the  same  part  within  the  same  model  can  be  replaced  by  a  simpler  geometry  in  one 
analysis  type,  but  not  within  another  analysis  type.  A  third  possible  area  for  future 
work  is  automatic  calculation  and  construction  of  a  simpler  part  geometry  for  the 
parts  capable  of  being  replaced  by  the  simpler  geometry. 
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